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Summary 

This  report  presents  research  work  perfonned  under  the  contract  No.  F30602-00-1-0515  with  the 
Air  Force  Research  Lab,  Rome,  NY.  This  report  consists  of  five  chapters,  each  containing 
introduction,  numbered  equations,  simulations,  figures,  conclusion,  and  references. 

In  Chaper  1,  we  consider  interference  suppression  and  multipath  mitigation  in  global  navigation 
satellite  systems  (GNSS).  In  particular,  a  self-coherence  anti-jamming  scheme  is  introduced  which  relies 
on  the  unique  structure  of  the  coarse/acquisition  (C/A)  code  of  the  satellite  signals.  Because  of  the 
repetition  of  the  C/A-code  within  each  navigation  symbol,  the  satellite  signals  exhibit  strong  self¬ 
coherence  between  chip-rate  samples  separated  by  integer  multiples  of  the  spreading  gain.  The  proposed 
scheme  utilizes  this  inherent  self-coherence  property  to  excise  interferers  that  have  different  temporal 
structures  from  that  of  the  satellite  signals.  Using  a  multi-antenna  navigation  receiver,  the  proposed 
approach  obtains  the  optimal  set  of  beamforming  coefficients  by  maximizing  the  cross-correlation 
between  the  output  signal  and  a  reference  signal,  which  is  generated  from  the  received  data.  It  is 
demonstrated  that  the  proposed  scheme  can  provide  high  gains  towards  all  satellite  signals  in  the  field  of 
view,  while  suppressing  strong  interferers.  By  imposing  constraints  on  the  beamformer,  the  proposed 
method  is  also  capable  of  mitigating  multipath  that  enters  the  receiver  from  or  near  the  horizon.  No 
knowledge  of  either  the  transmitted  navigation  symbols  or  the  satellite  positions  is  required. 

In  Chapter  2,  we  investigate  the  mitigation  of  frequency  modulated  (FM)  interference  in  GPS 
receivers.  In  difference  to  commonly  assumed  wideband  and  narrowband  interferers,  the  FM  interferers 
are  wideband,  but  instantaneously  narrowband,  and  as  such,  have  clear  time-frequency  (t-f)  signatures 
that  are  distinct  from  the  GPS  C/A  spread  spectrum  code.  In  the  proposed  technique,  the  estimate  of  the 
FM  interference  instantaneous  frequency  (IF)  and  the  interference  spatial  signature  are  used  to  construct 
the  spatio-temporal  interference  subspace.  The  IF  estimates  can  be  provided  using  existing  effective 
linear  or  bilinear  t-f  methods.  The  undesired  signal  arrival  is  suppressed  by  projecting  the  input  data  on 
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the  interference  orthogonal  subspace.  With  a  multi-sensor  receiver,  the  distinctions  in  both  the  spatial  and 
time-frequency  signatures  of  signal  arrivals  allow  effective  interference  suppressions.  This  chapter 
considers  the  deterministic  nature  of  the  signal  model  and  utilizes  the  known  underlying  structure  of  the 
GPS  C/A  code.  We  derive  the  receiver  SINR  under  exact  and  perturbed  IF  values.  The  effect  of  IF 
estimation  errors  on  both  Pseudo-range  measurements  and  navigation  data  recovery  is  analyzed. 
Simulation  results  comparing  the  receiver  performances  under  IF  errors  in  single  and  multi-antenna  GPS 
receivers  are  provided. 

In  Chapter  3,  we  continue  to  address  the  problem  of  suppressing  the  nonstationary  interference. 
Combined  spatial  and  time-frequency  signatures  of  signal  arrivals  at  a  multi-sensor  array  are  used  for 
nonstationary  interference  suppression  in  direct-sequence  spread-spectrum  (DS/SS)  communications. 
With  random  PN  spreading  code  and  deterministic  nonstationary  interferers,  the  use  of  antenna  arrays 
offers  increased  DS/SS  signal  dimensionality  relative  to  the  interferers.  Interference  mitigation  through 
spatio-temporal  subspace  projection  technique  leads  to  reduced  DS/SS  signal  distortion  and  improved 
performance  over  the  case  of  a  single  antenna  receiver.  The  angular  separation  between  the  interference 
and  desired  signals  is  shown  to  play  a  fundamental  role  in  trading  off  the  contribution  of  the  spatial  and 
time-frequency  signatures  to  the  interference  mitigation  process.  The  expressions  of  the  receiver  SINR 
implementing  subspace  projections  are  derived  and  numerical  results  are  provided. 

In  Chapter  4,  we  study  the  performance  of  the  delay  lock  loops  (DLL)  in  GPS  receivers  in  the 
presence  of  impulsive  noise.  The  use  of  GPS  has  broadened  to  include  mounting  on  or  inside  manned  or 
autonomous  vehicles  which  makes  it  subject  to  interference  generated  from  motor  emissions.  Many 
sources  of  interference  are  typically  modeled  as  impulsive  noise  whose  characteristics  may  vary  in  terms 
of  power,  pulse  width,  and  pulse  occurrences.  In  this  chapter,  we  examine  the  effect  of  impulsive  noise  on 
GPS  DLL.  We  consider  the  DLL  for  the  GPS  Coarse  Acquisition  code  (C/A),  which  is  used  in  civilian 
applications,  but  also  needed  in  military  GPS  receivers  to  perform  signal  acquisition  and  tracking.  We 
focus  on  the  statistics  of  the  noise  components  of  the  early,  late,  punctual  correlators,  which  contribute  to 
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the  discriminator  error.  The  discriminator  noise  components  are  produced  from  the  correlation  between 
the  impulsive  noise  and  the  early,  late  and  punctual  reference  C/A  code.  Due  to  long  time  averaging,  these 
components  assume  Gaussian  distributions.  The  discriminator  error  variance  is  derived,  incorporating  the 
front-end  precorrelation  filter.  It  is  shown  that  the  synchronization  error  variance  is  significantly  affected 
by  the  power  of  the  received  impulsive  noise,  the  precorrelation  filter,  and  the  sample  rate. 

Finally,  an  anti-jam  GPS  receiver  which  suppresses  interference  by  projecting  the  received  signal  on 
the  noise  subspace  obtained  via  subspace  tracking  is  proposed  in  Chapter  5.  The  resulting  interference- 
free  signal  is  then  processed  by  a  beamformer,  whose  weight  vector  is  obtained  by  maximizing  the  signal- 
to-noise  ratio  at  the  beamformer  output.  It  is  shown  that  the  proposed  receiver  can  effectively  eliminate 
interfere  and  enhance  the  GPS  signals  at  the  receiver  output. 
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Chapter  1 

A  Novel  Interference  Suppression  Scheme  for  Global 
Navigation  Satellite  Systems  Using  Antenna  Array 

1.  Introduction 

Satellite  navigation  is  a  tool  to  determine  position,  velocity,  and  precise  time  world  wide.  A  navigation 
receiving  device  determines  its  three  dimensional  position  plus  time  by  measuring  the  signal  traveling 
time  from  the  satellite  to  the  receiver  (the  so  called  pseudorange  due  to  the  clock  offset  at  the  receiver) 
[1],  The  generic  name  of  the  satellite  navigation  systems  is  Global  Navigation  Satellite  System  (GNSS). 
Currently,  there  are  two  operative  navigation  systems,  one  is  Global  Positioning  System  (GPS)  of  the 
United  States  and  the  other  one  is  Global  Navigation  Satellite  System  (GLONASS)  of  Russia  [2], 

Both  GPS  and  GLONASS  employ  direct  sequence  spread  spectrum  (DS  SS)  signaling.  The  GPS 
satellites  transmit  signals  at  two  L  band  frequencies:  L\  =  1. 57542GHz  and  L2  =  1. 2276GHz.  Each 
satellite  broadcasts  two  different  pseudorandom  (PRN)  codes,  a  coarse/acquisition  (C/A)  code  and  a 
precision  (P)  code,  using  code  division  multiple  access  (CDMA)  technique.  TheLl  carrier  transmits  both 
the  C/A  code  and  the  P  code,  whereas  the  L2  carrier  only  transmits  the  P  code  [1],  The  GPS  C/A  code  is  a 
Gold  code  with  a  chip  rate  of  1.023  Mchips/sec  (or  code  period  1023)  and  repeats  every  millisecond,  and 
the  P  code,  usually  encrypted  for  military  use,  has  a  chip  rate  10.23  Mchips/sec  and  repeats  about  every 
week.  Similar  to  GPS,  GLONASS  also  has  two  DS  SS  components.  However,  frequency  division 
multiple  access  (FDMA)  technique  is  used  in  GLONASS,  where  each  satellite  transmits  on  a  different 
center  frequency.  The  C/A  code  of  GLONASS  has  a  length  of  5 1 1  chips  at  a  chip  rate  of  5 1 1  kHz,  and  it 
is  the  same  for  every  GLONASS  satellite.  The  C/A  code  repeats  10  times  within  each  navigation  symbol 
which  has  a  rate  of  100  bps.  Another  component  of  GLONASS  has  10  times  the  chip  rate  (5.11  MHz)  of 
the  C/A  code  and  uses  a  longer  PRN  code.  In  this  chapter,  we  consider  only  the  signals  induced  by  the 
C/A  code. 
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Despite  the  ever  increasing  civilian  applications,  the  main  drawback  of  the  satellite  navigation 
systems  remains  to  be  its  high  sensibility  to  interference  and  multipath  [3],  [4],  [5],  [6],  [7],  [8],  which  are 
the  two  main  sources  of  errors  in  range  and  position  estimations.  The  effect  of  interference  on  the  GNSS 
receiver  is  to  reduce  the  signal  to  noise  ratio  (SNR)  of  the  navigation  signal  such  that  the  receiver  is 
unable  to  obtain  measurements  from  the  satellite,  thus  losing  its  ability  to  navigate  [1],  Jammers  reported 
to  impact  the  GPS  receivers  are  wideband  noise,  CW,  pulsed  noise,  pulsed  CW,  swept  tone  (chirped), 
frequency  hopping,  and  spoofers.  Each  type  of  jammers  has  advantages  and  drawbacks  in  terms  of 
complexity  and  effectiveness  [9].  The  spread  spectrum  (SS)  scheme,  which  underlines  the  GNSS  signal 
structure,  provides  a  certain  degree  of  protection  against  interference.  However,  when  the  interferer  power 
becomes  much  stronger  than  the  signal  power,  the  spreading  gain  alone  is  insufficient  to  yield  any 
meaningful  information.  For  example,  for  the  GPS  C/A  signal,  the  receiver  is  vulnerable  to  strong 
interferers  whose  power  exceeds  the  approximately  30  dB  gain  (101og10 1023  «  30  dB)  offered  via  the 

spreading/despreading  process.  It  is  thus  desirable  that  the  GNSS  receivers  operate  efficiently  in  the 
presence  of  strong  interference,  whether  it  is  intentional  or  unintentional. 

Interference  suppression  in  SS  communication  systems  has  been  an  active  research  topic  for  many 
years  and  a  number  of  techniques  have  been  developed  (see,  e.g.,  [4],  [10],  [11],  [12],  and  references 
therein).  In  satellite  navigation,  interference  can  be  combated  in  the  time,  space,  or  frequency  domain,  or 
in  a  domain  of  joint  variables,  e.g.,  time  frequency  [13],  [14]  or  space  time  [15],  [16].  Multiple  antenna 
receivers  allow  the  implementations  of  spatial  nulling  and  beamsteering  based  on  adaptive  beamforming 
and  high  resolution  direction  finding  methods.  These  methods  are  considered  to  be  effective  tools  for  anti 
jam  GPS  [9], 

Conventional  antenna  arrays,  which  are  only  based  on  spatial  processing,  are  among  the  simplest,  and 
yet  effective,  techniques  for  narrowband  interference  suppression.  Such  techniques,  however,  is 
inadequate  for  broadband  jammers  (such  as  spoofer)  cancellation  or  in  the  presence  of  multipath.  In  these 
cases,  the  temporal  degree  of  freedom  is  required.  Space  time  processing  provides  the  receiver  with 
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spatial  and  temporal  selectivity.  The  spatial  selectivity  allows  the  discrimination  between  the  navigation 
and  interference  signals  based  on  their  respective  Direction  of  Arrivals  (DOAs)  [17],  [18],  [19],  [20].  The 
temporal  selectivity  is  used  for  broadband  interference  and  multipath  cancellation.  Generally,  the  criteria 
for  determining  the  optimal  array  weights  include  maximum  signal  to  interference  plus  noise  ratio 
(MSINR),  minimum  mean  square  error  (MMSE),  and  minimum  output  power  (MOP)  [16].  The  MSINR 
approach  seeks  the  array  weight  vector  by  maximizing  the  receiver  output  SINR.  The  MMSE  method 
chooses  the  weight  vector  such  that  the  mean  square  difference  between  the  array  output  and  the  desired 
temporal  signal  is  minimized.  Since  the  navigation  signal  power  is  well  below  the  noise  floor  at  the 
receiver,  minimizing  the  output  power  while  attempting  to  preserve  the  navigation  signal  is  the  goal  of  the 
MOP  based  scheme.  While  these  methods  are  widely  used  in  interference  suppression  in  satellite 
navigation  systems,  one  obvious  drawback  is  that  they  all  require  some  kind  of  a  priori  knowledge  of  the 
problem  parameter  values.  For  example,  satellite  locations  are  needed  in  order  to  calculate  the  signal 
power  for  the  MSINR  and  MMSE  methods.  In  [21],  spatial  and  temporal  processing  techniques  are 
applied  to  remove  GPS  like  broadband  jammers  and  recover  the  navigation  information.  The  assumptions 
made  in  [21]  are  that  the  chip  and  bit  synchronizations  are  achieved,  implying  that  pseudorange 
measurements  are  obtained.  However,  the  assumption  of  satellite  positions  or  acquisitions  is  difficult  to 
enforce  under  persistent  jamming,  or  during  the  initial  satellite  searching  stage  when  any  synchronization 
is  yet  to  be  established. 

In  addition  to  interference,  GPS  pseudorange  and  carrier  phase  measurements  also  suffer  from  a 
variety  of  systematic  biases,  including  satellite  orbit  prediction  error  and  clock  drift,  ionospheric  and 
tropospheric  delay,  GPS  receiver  clock  offset,  and  signal  multipath  [22].  The  satellite  orbit  and  timing, 
ionospheric,  and  troposhperic  errors  can  be  removed  by  differencing  techniques  or  significantly  reduced 
by  modeling  [1].  The  receiver  clock  offset  can  also  be  corrected  by  differencing  but  is  often  solved  for  as 
an  unknown  in  the  position  solution.  Multipath,  on  the  other  hand,  is  normally  uncorrelated  between 
antenna  locations.  As  a  result,  differencing  will  not  cancel  the  errors  caused  by  multipath.  Also,  modeling 
multipath  for  each  antenna  location  is  difficult  and  impractical  [23],  To  combat  signal  multipath,  many 
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techniques  have  been  proposed.  Among  them,  narrow  correlator  is  one  of  the  most  widely  used 
approaches  that  improve  the  C/A  code  tracking  performance  by  reducing  the  space  between  the  early  and 
late  correlator  [24].  Other  multipath  mitigation  techniques  include  multipath  elimination  delay  lock  loop 
(MEDLL)  [25]  and  multipath  estimation  technology  (MET)  [22],  etc. 

This  chapter  proposes  a  new  interference  suppression  technique  for  GNSS  using  spatial  processing, 
but  incorporating  the  known  temporal  structure  of  the  C/A  signal.  A  careful  examination  of  the  existing 
interference  cancellation  techniques  reveals  that,  though  efficient  in  most  situations,  they  do  not  fully  take 
advantage  of  the  unique  C/A  signal  structure,  namely  the  replication  of  the  C/A  code.  Due  to  the 
repetition  of  the  spreading  code,  the  GNSS  C/A  signal  exhibits  strong  self  coherence  between  chip 
samples  that  are  separated  by  integer  multiples  of  the  spreading  gain.  Utilizing  this  feature,  an  anti 
jamming  technique  is  developed  to  suppress  a  large  class  of  narrowband  and  broadband  interferers.  It  also 
has  the  capability  of  mitigating  multipath,  resulting  in  improved  accuracy  in  pseudorange  measurements. 
The  proposed  technique  allows  the  civilian  C/A  code  tracking  and  acquisition  operations  in  the  presence 
of  strong  interference,  specifically  at  “cold  start”,  where  there  is  no  prior  information  on  satellite  angular 
positions  or  ranges  [26].  In  military  applications,  the  encrypted  P  code  is  used  instead  of  the  C/A  code. 
However,  due  to  the  short  duration  of  the  P  code  chip  (10  times  shorter  than  the  C/A  code  chip),  the 
synchronization  in  P  code  is  usually  difficult  to  achieve  using  the  early  late  correlator,  and  assistance 
from  the  C/A  code  is  needed  [24].  With  the  receiver  introduced  in  this  chapter,  initial  synchronization  in 
the  P  code  can  be  established  by  first  processing  the  interference  suppressed  C/A  code.  In  essence,  the  self 
coherence  based  anti  jamming  approach  is  a  blind  technique,  which  does  not  require  the  knowledge  of  the 
navigation  data  or  satellite  locations  to  perform  interference  suppression.  This  makes  it  most  applicable  in 
the  initial  satellite  searching  phase  when  such  information  is  unavailable,  or  in  a  prolonged  jamming 
environment  where  the  formerly  obtained  satellite  positions  are  no  longer  reliable. 

The  rest  of  the  chapter  is  organized  as  follows.  In  Section  2,  we  briefly  review  the  concept  of  spectral 
self  coherence  restoral.  An  antijamming  GNSS  scheme  is  presented  in  Section  3  and  two  receivers  based 
on  such  scheme  are  developed.  In  Section  4,  we  discuss  various  important  issues  related  to  the 
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performance  of  the  proposed  GPS  receiver.  Numerical  results  are  presented  in  Section  5  to  demonstrate 
the  performance  of  the  proposed  receiver.  Finally,  Section  6  concludes  this  chapter. 

2.  Overview  of  SCORE  Algorithms 

The  proposed  anti  jamming  GNSS  technique  builds  on  the  basic  concept  of  the  self  coherence  restoral 
algorithm  proposed  in  [27].  A  signal  s(t)  is  referred  to  as  spectrally  self  coherent  at  frequency  separation  fi 
if  the  correlation  between  the  signal  and  its  frequency  shifted  version  is  nonzero  for  some  lag  x,  i.e.,  if 
[27],  [28] 


where  (■)*  denotes  the  complex  conjugate  and  (i  'f  represents  the  infinite  time  averaging 

operation,  (r)  is  the  self-coherence  function  and  Rss( 0)  and  /Tf  ’(r)  represent  the  average  power 
and  cyclic  autocorrelation  function  of  s(t),  respectively.  For  an  M-element  vector  waveform  x(t),  the 
cyclic  autocorrelation  matrix  Rf'(r)  is  defined  as 

(x(t)xH(t-T)e-J2^)x  (1.2) 

where  (• )H  stands  for  the  complex  conjugate  transpose.  Complex  wide  sense  cyclostationary  waveforms 
exhibit  spectral  self  coherence  at  discrete  multiples  of  the  time  periodicities  of  the  waveform  statistics 
[27].  The  signal  waveforms  that  possess  the  self  coherence  feature  include  most  communication  signals, 
such  as  PCM  signals  and  BPSK  signals  [28]. 

The  spectral  self  coherence  restoral  (SCORE)  beamforming  techniques  have  been  shown  to  blindly 
extract  the  desired  signal  in  the  presence  of  unknown  noise  and  interference  [27].  The  SCORE  algorithms 
seek  the  beamformer  weight  vector  that  maximizes  a  measure  of  the  cyclic  feature  of  the  beamformer 
output.  For  example,  in  the  presence  of  interference,  the  received  signal  is  given  by  x  ( t )  =  as(  t)  +  v  ( t ), 
where  a  is  the  signal  amplitude  and  v  ( t )  is  the  interference,  which  is  assumed  to  be  independent  of  s(t).  If 
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s(t )  is  spectrally  self  coherent  at  a  frequency  shift  /?,  then  the  cyclic  autocorrelation  of  x(t)  can  be 
expressed  as  [27] 

=  |a|2  +  =  \4  (1.3) 

Equation  (1.3)  shows  that  the  shift  in  frequency  completely  decorrelates  the  interference  component  in 
x(t),  given  that  v{t)  is  not  spectrally  self-coherent  at  the  frequency  separation  fi. 

There  are  several  different  versions  of  the  SCORE  algorithm,  of  which  the  least-squares  (LS)  SCORE 
is  the  simplest.  The  LS-SCORE  algorithm  determines  the  array  weight  vector  by  minimizing  the 
difference  between  the  array  output  and  a  reference  signal,  which  is  obtained  by  processing  the  delayed 
and  frequency-shifted  version  of  the  received  signal.  Other  SCORE  algorithms  include  the  cross-SCORE 
algorithm,  which  determines  the  beamformer  by  strengthening  the  cross-correlation  between  the  output  of 
the  array  and  a  reference  signal,  and  the  auto-SCORE  algorithm,  which  maximizes  the  spectral  self¬ 
coherence  strength  at  the  output  of  a  linear  combiner  [27].  The  self-coherence  anti-jamming  scheme 
proposed  in  this  chapter  is  based  on  the  cross-SCORE  algorithm. 

3.  Proposed  Anti-Jamming  GNSS  Scheme 

Before  introducing  the  proposed  anti-jamming  receiver,  we  first  examine  the  temporal  structure  of  the 
navigation  signal,  as  the  receiver  is  developed  by  exploiting  the  repetitive  feature  of  the  C/A  signal. 
Figure  1  depicts  the  structure  of  the  received  noise-free  navigation  signal,  where  the  BPSK  modulated 
navigation  symbols  (simply  referred  to  as  “symbol”  thereafter)  are  spread  by  a  PRN  code  with  spreading 
gain  of  P  (P  =  1023  for  GPS  and  P  =  511  for  GLONASS)  and  chip-rate  sampled.  The  code  sequence 
(denoted  as  “spreading  block”  in  Figure  1)  is  repeated  L  times  ( L  =  20  for  GPS  and  L  =  10  for  GLONASS) 
within  each  symbol.  Two  blocks  of  data  are  formed  at  the  receiver:  a  data  block,  which  spans  N 
consecutive  samples,  and  a  reference  block  with  the  same  number  of  samples  as  the  date  block.  The 
distance  between  the  respective  samples  in  the  data  and  reference  blocks  is  set  equal  to  jP  chips,  where 
1  <  j  <L  .  Obviously,  due  to  the  repetition  of  the  spreading  code,  the  nth  sample  in  the  data  block  will 
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have  the  same  value  as  the  corresponding  «th  sample  in  the  reference  block,  providing  that  the  two 
samples  belong  to  the  same  symbol. 

From  the  temporal  structure  of  the  C/A  signal,  we  observe  the  inherent  self-coherence  between 
samples  in  the  data  block  and  the  reference  block,  due  to  the  repetition  of  the  spreading  code.  Based  on 
this  observation,  a  novel  anti-jamming  technique  is  developed  in  this  chapter,  and  Figure  2  shows  a  block 
diagram  of  the  proposed  receiver  with  this  technique.  In  Figure  2,  an  M-e lement  array  is  deployed.  There 
are  two  beamformers  in  the  receiver:  a  main  beamformer  w,  processing  samples  in  the  data  block,  and  an 
auxiliary  beamformer  f,  handling  data  from  the  reference  block.  An  error  signal  e(t)  is  formed  as  the 
difference  between  the  beamformer  output  z(t)  and  a  reference  signal  d(t),  which  is  the  output  of  f  .  For 
the  proposed  scheme,  the  weight  vectors  w  and  f  are  updated  according  to  the  cross-SCORE  algorithm. 

The  signal  reaching  the  GNSS  receiver  is  the  aggregate  of  the  satellite  navigation  signals,  their  respective 
multipaths,  additive  white  Gaussian  noise  (AWGN),  and  broadband/narrowband  interferers.  Thus,  after 
carrier  synchronization,  the  signal  received  at  the  receiver  can  be  expressed  as 

x(«)  =  YjS, (nT*  "  T \ )  +  Z Bkuk 0)d *  +  v(» 

“=0  Q  ^  K  (L4) 

=  (nTs  - z"o +  Z5? (nTs  ~ Tq )*/*'  +  Z Bkuk 0)d *  +  v(») 

q= 1  k=  1 

where  Ts  is  the  Nyquist  sampling  interval,  Q  is  the  number  of  multipath  components,  with  subscript  0 
designated  to  the  direct-path  signal.  Due  to  the  weak  cross-correlation  of  the  C/A-codes,  only  one  satellite 
is  considered  in  Equation  (1.4).  In  the  above  equation,  sq(n),  Tq  ,  and  (f)q  are  the  signal  sample,  time-delay, 

and  phase-shift  of  the  ^th  multipath  component,  respectively,  K  is  the  number  of  interferers,  uk(n)  is  the 
waveform  of  the  kth  interferer  with  amplitude  Bk.  The  vectors  aq  and  d*  are,  respectively,  M  x  1  spatial 
signatures  of  the  <yth  satellite  multipath  and  the  /cth  interferer,  and  \(n)  consists  of  noise  samples.  Let 

s(«)U  s0(nTs  -  r0  )a(leMl  denote  the  data  vector  across  the  array  due  to  the  direct-path  signal.  Then, 
Equation  (1.4)  can  be  rewritten  as 

x(»)  =  s(»)  +  s  (»)  +  u(»)  +  v(»)  (1.5) 
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Q  K 

where  s (n)  □  (nTs  -rq)aqeJ*q  and  u(n)  □  X  Bkuk (n)dk .  Assuming  a  direct-path  to  the  satellite  at 

q= 1  k= 1 

direction  ju  and  a  uniform  linear  array,  we  can  express  vector  a0  in  the  specific  format  of  a  steering  vector 
as 

a0  =  a(9)  □  [l,  eJ2xf‘T,  •••,  (1.6) 

where  fc  is  the  carrier  frequency,  v  =  A  /  c  sin  9  is  the  interelement  path  delay  of  the  source  in  the 
direction  of  6,  c  is  the  propagation  speed  of  the  waveform,  and  A  is  the  sensor  spacing.  According  to  the 
formulation  of  the  data  and  reference  blocks,  the  counterpart  of  \(n)  in  the  reference  block  within  the 
same  symbol  can  be  written  as 

x(«  -  jP)  =  YjS<,  (nTs  ~Tq~  jp)*/*q  +  X  Bkuk  0  -  jP) d*  +  v(«  -  jP) 

q= 0  k= 1 

=  XA/  in  T  -  vq  )aqeA  +  X  Bkuk  (n  -  jP)dk  +  \(n  -  jP )  (1.7) 

q= 0  k= 1 

=  s(n)  +  s (n)  +  u (n  -  jP )  +  \(n  -  jP ) 
where  we  have  assumed  that,  when  considered  within  the  same  symbol, 

Sq(nTs  ~Tq)  =  Sq(nTs  ~Tq  -  JP)’  <i  =  0>  "•»  Q  (L8) 

Compared  to  the  general  case  of  self-coherence,  there  is  no  frequency  difference  between  the  signal 
samples  in  the  data  and  reference  blocks,  i.e.,  the  frequency  shift  ft  =  0. 

From  Figure  2,  the  beamformer  output  and  the  reference  signal  are  given  by  z(n)  □  wH\(n)  and 
d(n)  D  f  Hx(n-jP) ,  respectively.  We  define  the  following  covariances: 

Rzd  □  E  | z(n)dH  (n)|  =  w H E  |x(«)x//  ( n  -  jP ) }  f  (1.9) 

R__  □  E {z(n)zH (n)}  =  wH E {x(ri)xH (n)}w  (1.10) 

Rdd  □  E{d(n)dH(n)}  =  fHE\x(n-jP)xH(n-jP)}f  (1.11) 

Under  the  assumption  that  the  navigation  signal,  interference,  and  noise  are  independent,  then 
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Rxx  □  E  (x(n)x77  («)}  =  E{x(n  -  jP)xH  (n  -  jP )}  =  R  +  R„  +  R 


(1-12) 


The  three  terms  on  the  right-hand  side  of  Equation  (1.12)  denote,  respectively,  the  covariance  matrices 
of  the  C/A  signal,  including  both  the  direct-path  and  multipath  signals,  interference,  and  noise: 

Rs  D  is  j[s(«)  +  s(«)][s(ft)  +  s(n)]^  j  ,  Ra  □  E  ju(«)u/y(n)j  ,  and  Rv  U  E {x(n)xH (n)}  .  If  the 

navigation  signal  is  the  only  data  component  which  correlates  with  its  delayed  version,  then  the  cross¬ 
correlation  matrix  between  the  corresponding  data  vectors  in  the  data  and  reference  blocks  simplifies  to 
[cf.  Equation  (1.3)] 

R^Dis{x(n)x>-yi>)}  =  Rs  (1.13) 


3.1.  Cross-SCORE  Algorithm  Based  Receiver 


We  first  consider  the  receiver  design  by  directly  applying  the  cross-SCORE  algorithm.  For  the  proposed 
GNSS  anti-jamming  scheme,  there  are  two  beamformers  w  and  f  to  be  determined.  With  d{n)  serving  as 

the  reference  signal,  we  define  e(n)U  z(n)-d(n )  as  the  difference  between  the  receiver’s  output  and 
the  reference  signal.  The  relationship  between  w  and  f  can  be  established  in  the  least-squares  (LS)  sense. 
For  a  fixed  beamformer  w,  the  LS  solution  of  f  is  given  by  fLS  =  R ,'rrx  ,  where 

r,v  =  E {x(n-  jP)zH («)}  =  R'f7/W .  Similarly,  if/is  fixed,  then  wLS  =  Rj'R^'f  ■ 

According  to  the  cross-SCORE  algorithm,  the  beamformers  w  and  f  are  obtained  by  maximizing  the 
cross-correlation  between  z(n)  and  d(n): 


C(w,f>8 


R 


zd 


w^R^f 


R:zRdd  [w"Rtrw][fffRv; 


(1.14) 


Substituting  f  and  w  in  the  above  equation  by  fLs  and  Wls,  respectively,  we  have 


wflR(P)R_1R(?)w 

f  \  _  XX  xx  xx 


w/7R  w 


(1.15) 
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C(wLS,f)  = 


(1.16) 


fH*J 

The  weight  vectors  w  and  f  that  maximize  C(w,  fLs)  and  C(wLS,  f),  respectively,  are  readily  shown  to  be 
the  eigenvectors  corresponding  to  the  largest  eigenvalues  of  the  generalized  eigenvalue  problems: 


R«w=zLrLp)r;Xp)w 

(1.17) 

=  irR^RjX^f 

(1.18) 

where  X  and  k  are  the  eigenvalues. 

It  is  observed  from  Equations  (1.17)  and  (1.18)  that,  for  the  proposed  receiver,  the  main  beamformer 
w,  which  generates  the  receiver  outputs,  is  equivalent  to  the  auxiliary  beamformer  f  that  provides  the 
reference  signal.  This  equivalence,  however,  is  not  surprising  because  of  the  unique  structure  of  the  C/A 
signal.  From  Equations  (1.4)  and  (1.7),  we  note  that  the  self-coherence  of  the  navigation  signal  is  due  to 
the  time  lag  between  the  two  samples  which  do  not  encounter  any  frequency  shift  after  the  frequency 
demodulation.  Therefore,  x(»)  and  \(n-jP)  have  the  same  correlation  function,  given  by  Equation  (1.11). 
For  the  general  case  of  self-coherence,  on  the  other  hand,  the  signal  auto-correlation  function  does  not 
necessarily  equal  to  the  auto-correlation  function  of  the  frequency-shifted,  time-lagged  version  of  the 
original  signal.  As  a  result,  the  cross-SCORE  algorithm  will  not  produce  two  identical  beamformers  [27], 
as  it  does  for  the  proposed  receiver. 

Simulations  presented  in  Section  V  show  that  the  cross-SCORE  based  receiver  performs  fairly  well  in 
a  jamming  environment.  It  is  capable  of  producing  high  gains  for  satellites  currently  in  the  field  of  view, 
while  suppressing  strong  jammers.  Elowever,  such  a  receiver  provides  no  measures  against  multipath, 
which  is  one  of  the  dominant  error  sources  in  navigation.  To  address  the  multipath  issue,  we  modify  the 
receiver  design,  as  discussed  in  the  next  subsection. 
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3.2.  Modified  Cross-SCORE  Algorithm  Based  Receiver 

Since  the  GNSS  signal  multipath  shares  the  same  structure  as  the  direct-path  signal,  it  is  expected  that 
both  the  C/A  signal  and  the  undesired  multipath  components  will  appear  at  the  receiver’s  output 
undistorted. 

It  is  known  that  the  satellites  typically  lie  above  the  horizon,  whereas  the  multipath  is  often  generated 
from  local  scatters  near  the  horizon  [16].  To  equip  the  receiver  with  means  to  combat  multipath,  while 
maintaining  the  self-coherent  approach,  we  introduce  constraints  on  f  such  that  the  reference  signal  d(n ) 
does  not  contain  reflections  from  near  the  horizon.  To  do  so,  we  define  equally  spaced  directions 
d=\,  -,D,  covering  some  solid  angle  Q  near  the  horizon.  Let  BU  [b(yj)  •••  b(y£))]  be  the  M  xD 

matrix  consisting  of  steering  vectors  defined  as  in  Equation  (6).  To  mitigate  multipath  in  the  range  Q ,  we 
require  Bwf  =  0.  Then,  the  cost  function  in  Equation  (1.16)  is  rewritten  as 

fopt  =  argmax  ,  subject  to  B  f  =  0  (1.19) 

f  1  Kt 

where  Rvv  =  R'^'R  j  R'  n// .  The  solution  of  the  above  equation  is  obtained  as  follows.  Let  r  =  rank(B) 

<  min (M,D)  be  the  rank  of  the  B  matrix.  Performing  the  singular  value  decomposition  (SVD)  [29]  of  B 
yields 

H  h  TA  0 

U"B"V=  (1.20) 

0  0 

where  U  and  V  are  two  unitary  matrices  with  dimension MxM and D  x  D,  respectively,  and 

A  =  diag{<r1,  a2,  •••  <rr}  (1.21) 

where  ax  >  <J2  >  ■  •  ■  >  cr.  are  eigenvalues  of  B  arranged  in  a  decreasing  order.  Let  A  be  formed  from  the 

last  M-r  columns  of  U.  Thus,  A  spans  the  null  space  of  BH,  i.e.,  BHA  =  0.  Let  a  be  a  ( M-r )  x  1  vector  such 
that 

f  =  Aa  (1.22) 
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Using  vector  a,  the  constrained  maximization  problem  in  Equation  (1.19)  is  transformed  to  an 
unconstrained  one: 


a  =  arg  max 

a 


aHAHRxxAa 

aHAHRxxAa 


(1.23) 


The  above  generalized  eigendecomposition  problem  can  be  solved  using  Cholesky  decomposition. 
Particularly,  since  Rxx  is  positive  definite,  A^R^A  is  also  positive  definite.  Then  the  Cholesky 
decomposition  of  AHRXXA  is  AhRxxA  =  GGh,  where  G  is  a  ( M-r )  x  ( M-r )  matrix  with  full  rank  [29]  and, 
thus,  invertible.  Let  a  =  G  77p  .  Then 

a"A"R„A«  P  C;  :A»R,AG  |1  p"G  A^R.AG  "p 
a"A"RvrAa  p"G  GG"G  "p  p"p 


Accordingly,  the  maximization  problem  becomes 


HAHR  Aa  p"G  ‘A"Rv  AG  "p 

—77 — 77  =  max - - 

a  A  RvvAa  p  p77 


=  maxp"G  A"Rv  AG  "p 

p 


(1.25) 


under  the  standard  constraint  ||p||  =  1,  where  ||-||  is  the  vector  2-norm  [29].  Hence,  P  is  given  by  the 
eigenvector  associated  with  the  maximum  eigenvalue  of  G  'A^R^.AG  // .  And,  finally, 


fopt=AG-/7p  (1.26) 

The  beamformer  w  is  derived  as 

«v  =  =  R;:RiTA«  (U7> 

according  to  the  LS  relation  between  w  and  f. 

4.  Covariance  Matrix  Estimations 


In  practice,  the  covariance  matrices  Rvvand  R,  n  are  unknown  and  have  to  be  replaced  by  their  sample 
estimates.  Define  the  M  x  N  data  and  reference  matrices  as  D  [x(«),  •••,  x(n-(Ar-l))]  and 
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XVrcf  □  [x(n  -  jP),  •••,  \(n  -  (A  -1)  -  jP)\ ,  where  N  is  the  block  length  and  N  <P.  The  sample 

covariance  matrices  are  then  given  by 

(1-28) 

(1.29) 

And,  the  beamformers  w  and  f  are  calculated  correspondingly. 

It  is  noted  from  Equation  (1.16)  that  the  covariance  matrices  R tT  and  R'^1  determine  the 
performance  of  the  proposed  receiver.  In  practical  implementations,  the  data  and  reference  blocks  X  v  and 
Xwref  are  used  to  estimate  Rxr  and  R1/^ ,  and  subsequently  provide  the  weight  vector  w,  which  is  then 
applied  to  process  signal  samples  in  the  data  block.  The  key  assumption  made  for  the  proposed  GPS 
receiver  in  Section  III  is  that  both  the  data  and  reference  samples,  x(n)  and  x(«-  jP),  1  <j<  20,  belong  to 
the  same  navigation  symbol.  However,  since  the  data  samples  used  for  covariance  matrix  estimations  are 
selected  randomly,  and  interference  suppression  is  performed  prior  to  any  symbol  synchronization 
process,  there  is  no  guarantee  that  the  data  and  reference  samples  belong  to  the  same  symbol.  Questions 
arise  as  how  will  the  receiver  perform  when  the  above  assumption  fails,  i.e.,  the  data  and  reference 
samples  lie  in  two  adjacent  symbols? 

To  answer  the  above  question,  we  relax  the  condition  imposed  in  Section  3,  and  develop  the  general 
expression  of  the  covariance  matrices  between  the  data  and  reference  samples,  x(«)  and  x(n-jP).  Define 
the  following  events: 

A  j :  x(n  )  &x(n-jP)  are  within  the  same  symbol, 

A  2 :  x(n)  &  x{n-  jP)  are  in  two  adjacent  symbols, 

A  21  :  x(n)  &x(n-  jP)  are  in  two  symbols  with  the  same  sign, 

A  22 :  x(n)  &  x(n  -  jP )  are  in  two  symbols  with  different  signs 

With  random  selection  of  time  n,  and  using  the  repetitive  property  of  the  C/A-code,  it  is 
straightforward  to  show  that  the  corresponding  probabilities  of  the  above  events  are 
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Pr  {A 1 }  =  T  T’P  =  1  - Jr .  Pr  {A  2 }  =  Jr  ’  Pr  iA  21 }  =  ^ »  and  Pr  {A  22 }  =  ^7 .  respectively,  where  T 

=  20 P  is  the  total  number  of  samples  in  one  symbol.  The  exact  expression  of  the  cross-correlation 
function  Rf  can  be  written  in  terms  of  the  above  probabilities  and  conditional  expectations  as 

R^1  =  E  jx(n)x"  (n  -  jP)\A  1 }  Pr  {A  J  +  E  {x(«)x"  (n  -  yP)|^2}Pr{^2} 

=  E{x(n)xH(n-jP)\Al)?x{Al}  (1.31) 

+E  jx(«)x"  (n  -  jP)\A  21}  Pr  {A  21 }  +  E  jx(«)x"  (n  -  jP)\A  22  J  Pr  [A  22 } 

Since  E {x(n)xH (n- jP)\A^  =  E ^x(n)xH (n  -  jP)\A2^  =  Rs  and 


E  |x(n)xw  (n  -  jP)\A  22 1  =  -R5 ,  then 


Rif  =R.s(Pr{^1}  +  Pr{^21})-RsPr{^22}  = 


\  T  j 


R 


(1-32) 


Equation  (1.32)  shows  that  the  covariance  matrix  R(f  depends  on  the  distance  between  the  data  and 


reference  samples  jP.  The  maximum  value  of  R(f  is  achieved  when  j  =  1,  representing  the  closest 
possible  data  and  reference  blocks. 

In  practice,  however,  sample  estimates,  obtained  from  Equations  (1.27)  and  (1.28)  using  the  data  and 
reference  blocks  XiV  and  XWref,  replace  the  exact  values  in  Equation  (1.32).  It  can  be  readily  shown  that  if 
XiV  and  X  \'rei  are  jP  samples  apart,  1  <  j  <  20,  the  probability  of  the  two  blocks  belonging  to  the  same 


symbol  or,  equivalently,  in  two  adjacent  symbols  with  the  same  sign,  is 


JP  +  N 
2  T 


On  the  other  hand. 


the  probability  that  X/V  and  X.vrei  are  in  two  adjacent  symbols  with  opposite  signs  is 


jP~N 

IT 


.  Using  the 


above  probabilities,  the  expected  values  of  Rv>  and  R*f  are  derived  in  Appendix  A  as: 


R  -R  R  R 


(1.33) 
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(1.34) 


which  show  the  same  dependency  on  jP  as  in  Equation  (1.32)  and  that  Rvvand  R^)  are  unbiased 
estimates  of  Rrx  and  R1^ ,  respectively. 

The  above  covariance  matrix  estimations  use  only  one  data  block  and  its  replicated  reference  block. 
To  fully  take  advantage  of  the  repetitive  feature  of  the  C/A-code,  multiple  data/reference  blocks  can  be 
used  in  the  time-averaging.  Particularly,  using  G  data  and  reference  blocks,  Equations  (1.27)  and  (1.28), 
respectively,  become 

R^=^ZX#,(g)X"(g)/JV  (1.35) 

G  g=i 

=4Zx„(g)X"nt(g)/A(  (1.36) 

G  g=l 

In  the  case  when  one  of  the  data  blocks  (and,  respectively,  a  reference  block)  is  split  between  two 
adjacent  symbols  with  opposite  signs,  a  maximum  of  only  two  of  the  G  terms  in  the  above  equation  may 
suffer  from  symbol  transition,  whereas  the  rest  of  the  terms  will  be  coherently  combined.  Appendix  A  and 
B  derive  the  mean  and  variance  of  the  above  estimations,  showing  the  value  of  using  a  higher  value  of  G. 

5.  Numerical  Results 

In  this  section,  we  evaluate  the  performance  of  the  proposed  self-coherence  anti-jamming  receiver  using 
the  GPS  C/A  signals. 

A  uniform  linear  array  (ULA)  consisting  of  M  =  7  sensors  with  half-wavelength  spacing  is  used  in 
simulations  with  one  satellite  and  no  multipath.  We  set  M-  11  for  simulations  with  multiple  satellites  or 
multipath.  The  GPS  navigation  symbols  are  in  the  BPSK  format  and  spread  by  C/A-codes  (Gold  codes) 
with  processing  gain  of  P  =  1023.  We  select  the  first  satellite  C/A-code  for  concept  demonstration.  At  the 
receiver,  chip-rate  sampling  is  performed  and  N  =  800  samples  are  collected  in  both  the  data  and 
reference  blocks  for  covariance  matrix  estimations.  The  signal-to-noise  ratio  (SNR)  and  signal-to- 
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interference-plus-noise  ratio  (SINR)  are  defined,  respectively,  as  SNR  =  1 0  logl0 1  /  Pv  and 
SINR  =  1 0  logl0 1  / (NP  +Pv),  all  in  dB,  where  unit  signal  power  is  assumed  for  simplicity,  Pv  is  the  noise 
power,  and  P,  is  the  interference  power.  The  jammer-to-signal  ratio  (JSR)  is  defined  as  JSR  =  101og10  P. 

dB.  Interferers  used  in  the  simulations  are  generated  as  broadband  binary  signals  having  the  same  rates  as 
the  C/A-codes,  but  with  a  different  structure  than  that  of  the  C/A  signals. 

5.1.  Antenna  Beam  Pattern  without  Interference 

We  first  consider  the  scenarios  in  which  no  interference  presents  at  the  receiver.  SNR  is  -30  dB. 
Covariance  matrices  are  estimated  using  one  data  block  and  one  reference  block  taken  within  the  same 
symbol.  The  performance  of  the  cross-SCORE  based  receiver  is  shown  in  Figure  3,  where  the  antenna 
pattern  is  formed  towards  the  satellite  located  at  6  =  30° . 

We  recall  the  discussions  in  Section  4  which  suggest  that  better  performance  can  be  expected  when 
multiple  data  and  reference  blocks  are  used  to  estimate  the  covariance  matrices.  In  Appendix  B,  the 
variances  of  the  sample  estimates  of  are  calculated  and  it  shows  that  using  multiple  data/reference 

blocks  can  indeed  reduce  the  estimation  variance.  We  now  demonstrate  experimentally  the  effect  of 
multi-block  estimation  on  the  receiver  performance.  Particularly,  a  very  special  situation  is  created  where 
one  of  the  data  blocks  is  evenly  split  between  two  symbols  having  opposite  signs.  SNR  is  set  at  -40  dB.  If 
G  =  2  and  the  split  data  block  happens  to  be  the  second  one,  the  receiver  fails  to  provide  any  substantial 
gain  for  the  satellite  located  at  30°,  as  shown  in  Figure  4(a).  This  is  because  that  elements  in  the  time¬ 
averaging  of  given  by  Equation  (1.35)  cancel  each  other,  resulting  in  significantly  weak  cross¬ 
correlation  between  z(n)  and  d(n)  [cf.  Equation  (1.14)].  If,  on  the  other  hand,  G  >  2  data  and  reference 
blocks  are  involved  in  the  estimation,  the  split  of  one  block  will  not  have  such  a  dramatic  impact  on  the 
receiver  performance  as  in  the  G  =  2  case,  as  only  two  among  G  blocks  are  affected  due  to  the  split.  It  is 
clear  from  Figure  4(b)  that  a  beam  is  generated  towards  the  satellite  with  G  =  7  despite  the  split. 
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Generally,  to  avoid  the  performance  degradation,  it  is  recommended  that  odd  number  of  samples  (N) 
should  be  chosen  for  each  data/reference  block  and  2  <  G  <  D. 

It  is  known  that  in  satellite  navigation,  at  least  four  satellites  are  needed  simultaneously  in  the  field  of 
view  in  order  to  calculate  the  receiver’s  three-dimensional  position  and  time.  Since  the  proposed  receiver 
relies  on  the  special  structure  of  the  C/A  signals  to  suppress  interference  and  all  satellite  emitted  C/A 
signals  share  the  same  repetitive  feature,  it  is  expected  that  the  receiver  will  pass  the  signals  from  all 
satellites  with  high  gains.  In  the  simulation,  the  satellites  are  located  at  6X  =10° ,  02  =30°,  6/  =50°, 

and  0A  =  70° ,  with  SNR  =  -30  dB.  As  shown  in  Figure  5,  four  clear  beams  are  generated  towards  the  four 
satellites. 

A  point  worth  mentioning  is  that  the  receiver  presented  in  this  chapter  is  able  to  suppress  interference 
for  all  satellites  at  once.  In  GNSS,  since  the  satellite  spreading  codes  are  known  at  the  receiver,  it  may  be 
intuitive  to  consider  using  the  spreading  code  as  the  reference  instead  of  generating  one  from  the  received 
data.  Even  though  the  locally  generated  spreading  code  is  noise  and  interference  free,  it  cannot  serve  as 
the  reference  signal  in  the  proposed  receiver  because  1)  the  alignment  of  the  incoming  signal  and  the  local 
reference  code  may  not  be  established  during  the  interference  removal  stage,  hence  there  is  no  guarantee 
that  the  data  block  and  reference  blocks  are  separated  by  integer  multiples  of  P  chips;  2)  it  is  not  possible 
to  use  one  specific  satellite’s  spreading  code  as  the  reference  signal  to  remove  interference  for  all 
satellites.  Therefore,  interference  suppression  must  occur  in  a  serial  manner.  The  disadvantage  is  obvious 
as  compared  to  the  simultaneous  interference  removal  the  proposed  receiver  offers. 

5.2.  Interference  Suppression 

We  next  investigate  the  receiver’s  interference  suppression  capability  by  comparing  it  with  the  MMSE 
receiver  of  [16].  The  MMSE  receiver  determines  the  weight  vector  by  minimizing  the  mean  square 
difference  between  the  array  output  and  the  desired  signal.  The  latter  approach,  however,  requires  the 
knowledge  of  the  satellite  direction.  This  condition  is  eliminated  in  the  proposed  scheme.  If  the  jammers 
have  explicit  bearings,  we  can  generate  the  received  signals  according  to  Equation  (1.4),  but  replacing  the 
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spatial  signature  d/  by  the  respective  steering  vector  defined  in  Equation  (1.6).  The  direction  of  the 
satellite  is  20° ,  while  two  jammers  are  located  at  40°  and  60° .  The  weight  vector  in  the  MMSE  method 
is  obtained  by  using  the  exact  transmitted  navigation  signal.  Figure  6  clearly  shows  that  deep  nulls  are 
placed  at  the  jammer  locations,  whereas  high  gains  are  generated  towards  the  direction  of  the  satellite  in 
both  schemes.  The  advantage  of  the  proposed  receiver  is  that  neither  prior  synchronization  nor  known 
satellite  location  is  required. 

5.3.  Multipath  Effects 

The  purpose  of  the  simulations  performed  in  this  subsection  is  to  demonstrate  the  difference  between  the 
receiver  that  solely  relies  on  the  cross-SCORE  algorithm  and  the  receiver  with  additional  constraint  in  the 
presence  of  multipath.  As  discussed  in  Section  3.1,  the  cross-SCORE  based  receiver  is  unable  to  mitigate 
the  signal  multipath,  though  very  efficient  in  suppressing  interference,  as  shown  by  simulations  presented 
so  far.  This  drawback  has  motivated  a  constrained  receiver  design  and  resulted  in  the  modified  cross- 
SCORE  algorithm  based  receiver  in  Section  3.2. 

We  assume  that  multipath  reaches  the  receiver  from  the  15-degree  range  (Q  as  defined  in  Section  3.2) 
above  the  horizon.  We  divide  Q  into  seven  2-degree  spaced  angles  and  form  the  corresponding  matrix  B. 
The  power  of  the  multipath  component  is  one  fifth  of  the  direct-path  signal  power. 

We  first  consider  the  case  when  there  is  no  interference.  The  direct-path  signal  is  incident  on  the 
array  with  50°  angle.  One  multipath  component  (half-chip  relative  delay  and  half  the  direct-path  signal 
power)  arrives  from  the  9°  direction.  Using  the  cross-SCORE  based  receiver,  both  the  direct-path  signal 
and  the  multipath  component  receive  high  gain  at  the  receiver  output  [Figure  7(a)].  If,  instead,  the 
modified  cross-SCORE  based  receiver  is  employed,  the  multipath  contribution  is  significantly  reduced 
from  the  output  of  the  receiver,  which  is  evident  from  Figure  7(b). 

In  the  next  case,  a  jammer  enters  into  the  system  from  30° .  Figure  8  shows  how  the  two  different 
receivers  respond  in  this  environment.  We  note  from  Figure  8(a)  that  both  receivers  can  successfully 
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place  deep  null  at  the  jammer  location.  However,  the  two  receivers’  responses  to  multipath  are  just 
opposite.  While  the  cross- SCORE  based  receiver  generates  a  beam  towards  the  multipath  component,  the 
modified  cross-SCORE  algorithm  based  receiver  creates  a  null  at  the  same  direction. 

The  multipath  mitigation  performance  of  the  proposed  receiver  is  also  evaluated  by  feeding  the  output 
of  the  receiver  to  a  conventional  early-late  delay  lock  loop  (DLL)  [24].  We  consider  the  discriminator 
functions  of  the  receiver  outputs  without  multipath  mitigation  and  the  outputs  with  the  modified  cross- 
SCORE  algorithm.  We  compare  the  results  with  the  case  where  there  is  no  multipath.  The  early-late 
spacing  is  set  to  be  half  of  the  C/A  chip  interval.  The  simulation  results  are  depicted  in  Figure  8(b),  which 
clearly  shows  that,  without  any  multipath  mitigation  process,  the  zero-crossing  point  of  the  discriminator 
function  drifts  away  from  the  origin,  indicating  the  pseudorange  measurement  error  [30].  If,  on  the  other 
hand,  the  modified  cross-SCORE  receiver  is  used  first  to  mitigate  multipath  contributions,  the  zero¬ 
crossing  point  of  the  corresponding  discriminator  function  almost  overlaps  with  the  zero-crossing  point 
obtained  using  the  direct-path  only  signal,  suggesting  that  the  proposed  technique  can  significantly  reduce 
the  multipath  effect  on  pseudorange  measurement. 

These  simulations  prove  that  both  receivers  have  the  capability  of  canceling  strong  jammers. 
However,  for  multipath  mitigation,  only  the  modified  cross-SCORE  algorithm  based  receiver  can  reject 
multipath  coming  from  near  the  horizon. 

5.4.  Synchronization  Process 

In  satellite  navigation,  the  receiver  is  ultimately  evaluated  by  its  ability  to  provide  accurate  pseudorange 
measurements.  This  is  achieved  by  establishing  synchronization  between  the  receiver  and  the  satellite, 
which  is  decided  based  on  the  cross-correlation  between  the  beamformer  outputs  and  a  locally  generated 
spreading  sequence  [31].  When  the  phase  of  the  receiver  replica  code  matches  that  of  the  code  sequence 
emitted  from  the  satellite,  there  is  a  maximum  correlation.  The  high-gain  beams  towards  the  satellites 
provided  in  the  previous  examples  should  be  examined  in  the  context  of  their  effects  on  the  post¬ 
processing  pseudorange  calculations. 
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In  the  simulation,  the  satellite  is  located  at  20°  and  the  two  jammers  are  at  40  and  60  .  The  figure 


of  merit  is  the  cross-correlation  between  the  receiver  output  and  the  Gold  code  sequence: 


CD 


E\  zHc 


£{Z"Zj£jc"cj 


(1.37) 


where  c  denotes  the  P  x  1  receiver  Gold  code,  z  is  a  P  x  1  vector  with  elements  given  by  z(n)  =  w " x 

and  w  is  the  beamformer  coefficient  vector  discussed  in  Section  3.1.  The  normalized  cross-correlation 
with  the  respective  antenna  beam  pattern  for  SNR  =  -25  dB  and  JSR  =  30  dB  and  50  dB  are  shown  in 
Figures  9  and  1 0,  respectively.  Also  shown  in  these  figures  are  the  normalized  cross-correlations  obtained 
before  the  jammers  are  removed.  It  is  observed  from  Figure  9(b)  that  synchronization  can  be  achieved  in 
the  presence  of  interference  when  JSR  =  30  dB.  Figure  10  shows  that  the  proposed  receiver  can 
effectively  cancel  directional  jammers  and  achieve  synchronization  even  when  the  JSR  is  as  high  as  50 
dB  [Figure  10(c)].  Without  interference  suppression,  however,  synchronization  fails  as  shown  in  Figure 
10(b). 


5.5.  Circular  Array 

In  addition  to  the  uniform  linear  array,  we  also  implemented  the  proposed  receiver  with  a  uniform  circular 
array  (UCA),  whose  configuration  is  shown  in  Figure  1 1(a).  Let  (0,  (f>)  denote  the  elevation  angle  and  the 
azimuth  angle  of  the  satellite.  Then,  the  steering  vector  of  the  satellite  for  the  M-element  UCA  is  given  by 


a  (0,0)  = 


/2  nr  sin  0  cos  d) 
Q  X 


.2 nr  .  Q  f  ,  0  M- 1 

J - Sint/COS  <p-L7t - 

_  X  l  M 


(1.38) 


where  r  is  the  radius  of  the  circular  array  and  /.  is  the  wavelength.  The  steering  vector  of  the  jammer  has 
the  same  form  as  a(  9,  (j)  given  above.  In  the  simulation,  the  satellite  signal  reaches  the  array  from 
(10° ,  20° ),  whereas  a  jammer  is  located  at  (  60° ,  40° ).  The  beam  pattern  is  shown  in  Figure  1 1(b)  for  r  = 
X,  and  SNR  =  -30  dB  and  JSR  =  30  dB.  It  is  observed  from  Figure  11(b)  that  the  receiver  has  the  ability  to 
reject  jammers  from  arbitrary  directions. 
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6.  Conclusions 


In  this  chapter,  we  addressed  the  issue  of  interference  suppression  in  global  satellite  navigation  system. 
Specifically,  the  unique  structure  of  the  GNSS  C/A  signal  is  exploited.  Due  to  the  repetition  of  the  C/A- 
code  within  each  navigation  symbol,  strong  self-coherence  is  observed  between  chip-rate  sampled  signals. 
It  is  shown  that  the  use  of  this  self-coherence  feature  allows  the  development  of  an  anti-jamming  GNSS 
receiving  scheme  which  is  built  on  the  cross-SCORE  algorithm.  The  proposed  scheme  incoiporates 
multiple  data  and  reference  blocks,  separated  by  integer  multiples  of  the  spreading  code  length,  to 
generate  the  array  weight  vectors.  Its  performance  is  analyzed  in  view  of  its  dependency  on  the  number 
of  blocks  and  the  number  of  samples  in  each  block.  Two  receivers  are  constructed  based  on  the  proposed 
scheme.  One  directly  applies  the  cross-SCORE  algorithm  which  seeks  the  optimal  beamformers  by 
maximizing  the  cross-correlation  between  the  receiver  output  and  a  reference  signal,  derived  from  the 
receiver  signal.  The  other  one  applies  constraints  on  the  beamformer  such  that  it  can  also  reject  multipath 
arriving  from  and  near  the  horizon.  Numerical  results  have  shown  that  the  proposed  scheme  is  capable  of 
suppressing  strong  wide  class  of  narrowband  and  broadband  interferers  while  preserving  signals  and  no  a 
priori  knowledge  of  either  the  transmitted  signals  or  the  satellite  locations  is  required. 
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Appendix  A  Mean  Calculation 

To  simplify  the  derivation,  we  rewrite  the  received  signal  vector  as 

x(«)  =  s(n)a(B)  +  \(n)  (A.l) 

where  we  consider  only  the  direct-path  signal  with  an  explicit  direction  dof  the  satellite.  The  vector  v(n) 
contains  samples  of  interference  and  noise,  with  zero  mean  and  variance  cr“ .  Both  components  of  \(n) 
are  assumed  to  be  independent  of  the  GPS  signal.  Accordingly, 
Rs  =  E  {s(n)a(0)aH  (0)sH  (n)}  =  a(6)an  (6)  □  Ra,  where  it  is  assumed 

£{5(n>//(«)}  =  £{|5(«)|2J  =  l. 

From  Section  3,  the  estimates  of  the  covariance  matrices  Rvvand  R'j'1  are  obtained  using  the  data 
and  reference  blocks  Xv  and  XVref  as  Rt;c  =  XvX!j / N  and  R'^-  X  vX^rcr//V  ,  respectively.  Taking 
the  expected  value  R  vv  yields 

A,  =  e{rJ =Tgir{x(„_0x*(„_0} = r, +R, 

N  1=0 

/  p\ 

The  expected  value  of  Rvv  is  obtained  as  follows.  Define  the  following  events: 

B  t  :  X  v  &  X;Vrcf  are  within  the  same  symbol, 

B2 :  X  v  or  XVrcf  is  split  between  two  adjacent  symbols 
&  iVj  <  N  samples  are  in  the  first  symbol, 

B 2 :  the  entire  X  v  and  X  Vrcf  are  in  two  adjacent  symbols. 


(A.2) 


(A.3) 
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The  corresponding  probabilities  of  the  above  events  are:  Pr {B , }  =  1  -  ^  ^  ,  Pr =  and 

Pr{2?3}  =  — — — ,  respectively.  In  addition,  we  define  the  following  events  regarding  the  two  adjacent 


symbols: 


Cj :  the  two  adjacent  symbols  have  the  same  sign, 
C2 :  the  two  adjacent  symbols  have  opposite  signs. 


(A.4) 


The  expected  value  of  is  calculated  as 


R 


(P) 


I  N  — 

=  E  {r”}  =  E  {x(n  -i)x"(n  -  i  -  jP )}  =  R 


(P)  .n^)  ,d(^) 

xx\Bx  '  ‘tVxx|52  ' 


(A.5) 


where 


1  N 


RS,  D  -Z^{x(»-0x/f(n-/-yP)|i?1}Pr{51 
jy  i= i 


yP  +  A 


1- 

v  T  j 


R 


(A.6) 


Because  the  GPS  symbols  are  equi-probable,  then  the  occurrence  of  B,  implies 


1  N 

RS,  □  T7  X E  *("  -<>"(«-  /  ■ -  ;/>)|  B,  Pr  {B, }  =  0 

i= i 


(A.7) 


On  the  other  hand,  when  event  B2  occurs, 


RL,k  □ 


i  N 


-YjE{TL{n-i)TLH{n-i-jP)\B1}^{B2}  =  -K 


A, 

T 


(A.8) 


From  Equations  (A.5  -  A.8),  the  expected  value  of  R^  is  given  by 


R 


( p > 


1- 


JP+N 


R,+yR, 


\_jP} 

T  j 


R 


(A.9) 


which  is  exactly  the  one  shown  in  Equation  (1.32). 

When  using  G  data  and  reference  blocks,  the  estimate  of  the  covariance  matrix  R'jy1  is  given  by 


RS=^Zx„(s)xlr(g)/v 

G  g=l 


(A.  10) 
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*  /  p\  'X  /  p\ 

To  calculate  R  ti(-; ,  the  expected  value  of  R  vr(l ,  we  define  the  following  events: 

F j :  the  first  data  block  or  the  last  reference  block  is  split, 

Fz:  one  of  the  other  G  - 1  data/reference  blocks  is  split 

&  Gj<  G  data  blocks  (respectively,  Gj  - 1  reference  blocks)  are  in  one  symbol,  (A.  11) 
F 3 :  the  data  blocks  and  reference  blocks  are  within  the  same  symbol, 

F4  :  no  split  block  and  the  data  and  reference  blocks  are  in  two  adjacent  symbols. 


The  corresponding  probabilities  are  Prj/'j] 


N_ 

T 


Pr{F2} 


N_ 

T 


Pr{F3}  =  l 


N  +  GjP 
T 


and 


Pr  {F4}  =  —  — - .  We  maintain  that 


=^S£{xUg)Xlf(g)|F,}Pr{F,}/Af=  t  c 


N  °~lR,+yR, 


(A.  12) 


and 


T*(F) 

,XxxG|F. 


=  ^  £ E  {: X„  (g)X“  „  (g  )|  F2  J  Pr  {F, }  /  IV  =  tfAA  PA  R <  +  (C_M_  R  ( A.  1 3) 


Further, 


r(F) 

rVxvG|F3 


=  -fZ£{x»<«)XL,(«)|FJ}Pr{FJ)/lV  =  fl  - 

vr  g=\ 


Rs  (A- 14) 


and 


D(f) 

rVxxG|F4 


=  FT£{X»(g)X"-(«)lF4PrT- 

(J  g=i 


IN 


GjP -N  G -2 

- R  + 

2  T  G 


GjP-N 

IT 


R.s  (A.  15) 


Finally,  R^  is  given  by 


A(f> 

atG 

_  n  (P)  i  p  (P)  i|3(^)  i  p  (p) 

^xxGIF,  ~r  m*xxG\F2  jvxcG|F3  "r  jvxxG|F4 

f  iP) 

=  I-*-  R? 
l  T  j 

(A.  16) 

which 

is 

equivalent 

to  the  expected  value 

given  in  Equation 

(1.32). 
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Appendix  B  Variance  Calculation 


The  variance  of  R(v^'  is  given  by  [32] 


var[R«]  =  £{RLPIR™}-£2{] 


(A.  17) 


In  Appendix  A,  we  have  shown  that  the  covariance  estimates  are  unbiased.  In  what  follows,  we 


concentrate  on  evaluating  E  j 


R(nR(m 


Using  one  data  and  reference  block,  we  let  □  E  j  j-  Pr  {/?  j ,  where  events  Bl , 


/e[l,3],  are  defined  in  Equation  (A.3).  Then, 


_|_u»  +\p 

1|S,  ^  1| B2  ^  1 1|S3 


(A.  18) 


When  the  data  and  reference  blocks  are  within  the  same  symbol  or,  equivalently,  in  two  adjacent  symbols 
with  the  same  sign  (i.e.,  event  C] ),  we  have 


E {R^R^  Q }  =  jg  g  x(n  -  i)xH (n  -  i  -  jP)x(n  - 1  -  jP)xH  (n  -  /)|  Q  j  (A.  1 9) 


Substituting  x(/i)  from  Equation  (1.39)  and  after  some  straightforward  calculations,  we  have 


XX  I M 


E 1  R(P)RiP)H | q  j,  =  M(N  +  Gv ) R  +M(l  +  o-„)R  qr 

N  N 


(A.20) 


where  we  have  used  a 77 (0)a(0)  =  M  and  E  jv77(»)v(»)}  =  Mcr2  ■  When  event  Bl  occurs, 


%.r,.z±^jR 

V  1  J 


(A.21) 


Similarly, 


2  N  ~ 

= - R 


(A.22) 


and  in  case  of  event  B, ,  we  have 
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(A.23) 


*F  =  {jP  — R 

l|i>3  rji 


Finally,  the  variance  of  R^1  is 


var[R«]  =  £{R«R'f“}-£2{] 


M(A'+ct,;)r  |  M(l  +  CT;)n  _ri_7P'|  r! 


(A.24) 


For  G  data  and  reference  blocks  we  can  similarly  define,  using  the  events  in  Equation  (A.l  1), 


(A.25) 


where  □  ii iR^gR^7  ,  /e[l,4],  Following  the  same  procedure  we  adopted  in 


calculating  the  expected  value  of  ,  it  can  be  readily  shown  that 


TG|F, 


N(G-\f  |  N  ^ 

T  G2  T  \ 

(G-Y)N  (G-l)2  |  (G-1)AT~L 
2  T  G2  2  T 

i_y-QpV 


G(,/P-iV)(G-l)2  |  GIJP-N)  R 


(A.26) 


from  which  we  obtain 


E\  R;R 


,-2Z+2  Z  1  ^A)R  +m±aj)R 


(A.27) 


Finally,  the  variance  is 


var[R^]  =  £{R«R™}-£2{] 


=fi_2Z+2A'|  ZEZr  ,  M(i+g;)tt 

v  TGT  N  “  N 


;P\2  (A.28) 

1-^-  Rl 
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The  above  equation  clearly  shows  that  the  larger  the  number  of  data  and  reference  blocks  used  in  the 


time-averaging,  the  smaller  the  variance. 
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Figure  2.  Structure  of  the  proposed  anti-jamming  scheme. 
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Figure  3.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  with  SNR  =  -30  dB  and  one 
data  and  one  reference  block  taking  within  the  same  symbol. 

90  90 

1  1 


(a) 


(b) 


Figure  4.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  with  multiple  data  and 
reference  blocks  and  SNR  =  -40  dB.  (a)  G  =  2;  (b)  G  =  7. 
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Gain  in  dB 


Figure  5.  Beam  pattern  generated  by  the  cross-SCORE  based  receiver  with  four  satellites,  SNR  = 

30  dB,  and  (7=7. 


Figure  6.  Antenna  gains  of  the  proposed  scheme  and  the  MMSE  scheme  with  SINR  =  -33  dB,  JSR 

30  dB,  and  G  =  3. 
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(a)  (b) 

Figure  7.  In  the  presence  of  multipath  with  SNR  =  -30  dB  and  G  =  l.  (a)  Cross-SCORE  based 
receiver;  (b)  Modified  cross-SCORE  based  receiver. 


(a)  (b) 


Figure  8.  (a)  Comparison  between  the  cross-SCORE  based  receiver  and  the  modified  cross-SCORE 
based  receiver  SINR  =  -  33  dB,  JSR  =  30  dB,  and  G  =  7;  (b)  Comparison  of  the  discriminator 

functions  of  the  early-late  delay  lock  loop. 
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Figure  9.  Synchronization  of  the  cross-SCORE  based  receiver  with  SNR  =  -25  dB,  JSR  =  30  dB, 
and  (7  =  7.  (a)  Beam  pattern;  (b)  Normalized  cross-correlation  before  jammer  removal;  (c) 
Normalized  cross-correlation  after  jammer  removal. 
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Figure  11.  Performance  of  the  receiver  with  a  circular  array,  (a)  Circular  array  configuration;  (b) 

Beam  pattern. 
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Chapter  2 

Subspace  Array  Processing  for  the  Suppression  of  FM 
Jamming  in  GPS  Receivers 


1.  Introduction 

Signals  with  rapidly  time-varying  frequency  characteristics  can  be  the  preferred  option  at  different 
transmitters  or  generated  as  a  consequence  of  moving  targets  and  scatterers  of  fixed-frequency  emitters 
[12].  These  signals  are  difficult  to  analyze,  track,  and  remove  using  conventional  methods  assuming 
stationary  or  quasi-stationary  environments.  Recent  developments  in  time-frequency  signal 
representations  [9],  [13],  [25]  have  allowed  accurate  estimations  of  the  signal  time-frequency  signature, 
including  the  signal  instantaneous  frequency  and  instantaneous  bandwidth,  leading  to  effective 
suppression  of  undesired  nonstationary  signals  [1],  [5],  [15],  [27],  Previous  work  has  addressed  the 
general  case  of  direct  sequence  spread  spectrum  receivers.  In  this  chapter,  we  consider  the  specific  case  of 
suppressing  a  class  of  frequency-modulated  jammers  incident  on  a  Global  Position  System  (GPS)  receiver 
using  the  jammer  spatial  and  time-frequency  signatures  in  conjunction  with  methods  of  subspace 
projections. 

The  GPS  uses  a  direct  sequence  spread  spectrum  (DSSS)  signal  that  is  highly  susceptible  to 
interference  [16],  [24],  There  are  both  intentional  and  non- intentional  forms  of  interference.  The  primary 
area  of  concern  in  both  commercial  and  military  application  is  intentional  interferers.  Development  of 
techniques  for  protection  of  GPS  from  interference  and  jamming  is  an  area  of  active  research. 

GPS  DSSS  signals  have  some  degree  of  jamming  protection,  via  processing  gain,  built  in  to  the  signal 
structure  itself;  however,  due  to  the  fact  that  the  GPS  signal  originates  in  a  half-geosynchronous  orbit,  it  is 
relatively  weak  when  it  reaches  the  earth.  The  weak  signal  strength  of  the  GPS  signal  makes  it  easy  for  an 
intentional  jammer  to  overcome  the  inherent  jamming  protection  of  the  DSSS  signal. 
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Depending  on  the  GPS  receiver,  the  jammer-to-signal  ratio  (JSN)  of  greater  than  40  to  50dB  will 
prevent  the  GPS  receiver  from  being  able  to  obtain  a  position.  There  are  several  approaches  to  mitigating 
this  susceptibility  including  frequency-domain  techniques  [6],  [11],  time-domain  techniques  [17],  [21], 
and  adaptive  antennas  [14],  [23],  Frequency-  and  time-domain  techniques  are  only  effective  against 
partial  band  interference.  Further,  neither  of  these  techniques  is  capable  of  effectively  incorporating  the 
suddenly-  changing  or  evolutionary  rapidly  time-varying  nature  of  the  frequency  characteristics  of  the 
interference.  In  both  techniques,  there  is  a  lack  of  intelligence  about  interference  behavior  in  the  joint 
time-frequency  (t-f)  domain,  rendering  them  limited  in  results  and  applicability.  For  the  time-varying 
interference  depicted  in  Figure  1,  frequency-domain  methods  remove  the  frequency  band  Af  and  ignore 
the  fact  that  only  few  frequency  bins  are  contaminated  by  the  interference  at  a  given  time.  Dually,  time- 
domain  excision  techniques,  through  gating  or  clipping  the  interference  over  AT ,  do  not  account  for  the 
cases  where  only  few  time  samples  are  contaminated  by  the  interference  for  a  given  frequency.  Applying 
either  method  will  indeed  eliminate  the  interference  but  at  the  cost  of  unnecessarily  reducing  the  desired 
signal  energy. 

Adaptive  excision  methods  might  be  able  to  track  and  remove  the  nonstationary  interference,  but 
would  fail  if  the  interference  is  highly  nonlinear  FM  or  linear  FM,  as  in  Figure  1,  with  high  sweep  rates. 
Further,  the  adaptive  filtering  length  or  block  transform  length  trades  off  the  temporal  and  the  spectral 
resolutions  of  the  interference.  Increasing  the  step  size  parameter  increases  the  filter  output  errors  at 
convergence,  and  causes  an  unstable  estimate  of  the  interference  waveform. 

The  above  example  clearly  demonstrates  that  nonstationary  interferers,  which  have  model  parameters 
that  rapidly  change  with  time,  are  particularly  troublesome  due  to  the  inability  of  single-domain 
mitigation  algorithms  to  adequately  ameliorate  their  effects.  In  this  challenging  situation,  and  others  like 
it,  joint  t-f  techniques  can  provide  significant  performance  gains,  since  the  instantaneous  frequency  (IF), 
the  instantaneous  bandwidth,  and  the  energy  measurement,  in  addition  to  myriad  other  parameters,  are 
available.  The  objective  is  then  to  estimate  the  t-f  signature  of  the  received  data  using  t-f  analysis, 
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attenuating,  or  removing,  the  received  signal  in  those  t-f  regions  that  contain  strong  interference.  This  is 
depicted  by  the  region  in  between  the  dashed  lines  in  Figure  1 . 

Recently,  several  techniques  based  on  linear  transforms  and  quadratic  distributions  have  been  devised 
for  FM  interference  excision  in  direct-sequence  spread-spectrum  (DSSS)  communication  systems  [1]  - 
[5],  [15],  [18],  [19],  [27],  [28].  These  techniques  assume  clear  jammer  time-frequency  signatures  and  rely 
on  the  distinct  differences  in  the  localization  properties  between  the  jammer  and  the  spread  spectrum 
signals.  An  important  class  of  these  t-f  based  methods  applies  subspace  projection  techniques  for 
interference  mitigation  [4],  [26],  [31],  [32].  In  essence,  the  jammer  instantaneous  frequency,  whether 
provided  by  the  time-frequency  distributions  or  any  other  IF  estimator,  is  used  to  define  the  temporal 
signature  of  the  interference,  with  one-dimensional  signal  space  per  interference  source.  This,  in  turn,  is 
used  to  construct  a  subspace  orthogonal  to  the  jammer.  The  respective  projection  matrix  is  used  to  excise 
the  jammer  power  in  the  incoming  signal  prior  to  correlation  with  the  receiver  pseudorandom  noise  (PN) 
sequence.  The  result  is  improved  receiver  signal-to-interference-plus-noise  ratio  (SINR)  and  reduced 
BERs.  Compared  with  the  subspace  projection  techniques  in  the  single-sensor  case,  the  use  of  multi¬ 
sensor  array  greatly  increases  the  dimension  of  the  available  signal  subspace.  It  allows  both  the 
distinctions  in  the  spatial  and  time-frequency  (temporal)  signatures  of  the  GPS  signals  from  those  of  the 
interferers  to  play  equal  roles  in  suppressing  the  jammer  with  a  minimum  distortion  of  the  desired  signal. 

In  this  chapter,  we  examine  the  applicability  of  multi-sensor  subspace  projection  techniques  for 
suppressing  nonstationary  jammers  in  GPS  receivers.  We  rely  on  IF  estimators,  specifically  the  TFDs,  to 
provide  accurate  estimates  of  the  jammer  signal  parameters.  With  any  employed  IF  estimator,  the  jammer 
parameters  are  deemed  to  be  perturbed.  Due  to  inaccuracies  in  IF,  the  GPS  receiver  anti-jamming 
performance  is  degraded,  lowering  the  receiver  SINR.  This  chapter  analyzes  the  multi-antenna  GPS 
receiver  performance  in  the  presence  of  zero-mean  identical  and  independent  Gaussian  IF  estimation 
errors.  The  single  antenna  receiver  case  is  derived  as  a  special  case  of  the  multi-antenna  receiver.  It  is 
shown  that  the  use  of  several  antennas  at  the  receivers  reduce  the  impact  of  IF  estimation  errors  on 
degrading  the  receiver  performance.  The  angular  location  of  each  GPS  satellite  is  assumed  known  and  is 
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used  to  design  an  appropriate  jammer  mitigation  technique.  In  addition,  accurate  estimates  of  the  jammer 
spatial  signatures  are  assumed,  and  stem  from  recently  developed  successful  direction  finding  techniques 
of  FM  signals  [7],  [8],  [30].  In  comparing  the  single  and  multi-antenna  cases,  the  chapter  shows  that  the 
use  of  antenna  arrays  effectively  improves  the  receiver  SINR  by  exploiting  the  difference  in  spatial 
signatures  as  well  as  the  t-f  signatures. 

This  chapter  is  organized  as  follows.  Section  2  gives  a  description  of  the  GPS  signal  structure.  It  also 
briefly  discusses  nonstationary  signal  parameter  estimations,  and  provides  key  references  on  the  subject. 
FM  jammer  suppression  using  subspace  projection  is  introduced  and  analyzed  in  Section  3  for  multi¬ 
antenna  GPS  receivers.  The  effect  of  inaccuracies  in  IF  estimation  on  receiver  performance  is  analyzed  in 
Section  4  for  the  single -jammer  case.  Appendix  A  presents  the  generalization  of  the  analysis  of  Section  4 
to  multiple  jammers.  All  jammers  are  modeled  as  frequency  modulated  signals  with  no  instantaneous 
bandwidth.  AM/FM  jammers  are  outside  the  scope  of  this  chapter  and  require  a  different  analytical 
approach. 

2.  Background 

2.1.  GPS  C/A  Signal  Structure 

GPS  employs  BPSK-modulated  DSSS  signals.  The  navigation  data  is  transmitted  at  a  symbol  rate  of  50 
bps.  It  is  spread  by  a  coarse  acquisition  (C/A)  code  and  a  precision  (P)  code.  The  C/A  code  is  a  Gold 
sequence  with  a  chip  rate  of  1.023  MHz  and  a  period  of  1023  chips,  i.e.  its  period  is  1  ms,  and  there  are 
20  periods  within  one  data  symbol.  The  P  code  is  a  PN  code  at  the  rate  of  10.23  MHz  and  with  a  period 
of  1  week.  These  two  spreading  codes  are  multiplexed  in  quadrature  phases[l 0,1 6,24],  Figure  2  shows 
the  signal  structure.  The  carrier  L 1  is  modulated  by  both  the  C/A  code  and  the  P  code  and  the  carrier  L2  is 
only  modulated  by  P  code.  It  is  commonly  assumed  that  the  C/A  code  and  the  P  code  are  perfectly 
separated.  The  peak  power  spectral  density  of  the  C/A  signal  exceeds  that  of  the  P  code  by  13dB.  The 
GPS  signal  is  typically  very  weak.  The  jammer-to-signal  ratio  (JSR)  is  often  larger  than  40  dB,  whereas 
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the  signal-to-noise  ratio  (SNR)  is  about  -14  to  -20dB.  Due  to  the  high  JSR,  the  jammer  depicts  a  clear 
signature  in  the  time-frequency  domain. 

We  mainly  address  the  problem  of  anti-jamming  for  the  C/A  code.  This  code  has  unique  properties 
due  to  the  fact  that  it  is  periodic  within  one  data  symbol  and  fixed  for  each  satellite  signal.  It  is  shown 
that  this  periodicity  has  a  considerable  effect  on  the  receiver  performance. 

2.2.  Instantaneous  Frequency  Estimation 

As  stated  previously,  there  are  several  linear  and  bilinear  t-f  methods  for  IF  estimation,  most  notably, 
those  based  on  Wigner-Ville  distribution  and  its  generalization  into  Cohen’s  class.  The  Wigner 
distribution  (WD)  (also  known  as  Wigner-Ville  distribution),  which  was  the  first  distribution  introduced 
in  the  context  of  quantum  mechanics,  has  paved  the  way  to  several  key  contributions  to  advances  in  the 
area  of  time-frequency  analysis  as  well  as  representations  of  signals  with  time-varying  characteristics. 
These  contributions  have  aimed  at  overcoming  the  drawbacks  of  the  WD  and  sought  new,  more  effective 
tools  for  nonstationary  signal  analysis,  synthesis,  and  processing. 

Cohen  [13]  provided  a  consistent  set  of  definitions  for  a  desirable  class  of  time-frequency 
distributions  (TFDs),  often  referred  as  Cohen’s  class.  Cohen’s  class  of  time-frequency  (t,co)  distributions 
for  the  signal  x(t)  may  be  presented  in  different  forms,  including 

co  oo 

pit, co; (/>)  =  J  |  <f>(t -u,t)x(u  +  t/2)x(u -t/2)  e~iC0Tdu  dr  (2.1) 

—oo  —oo 

Different  distributions  are  obtained  by  selecting  different  kernels,  (j)(t,T)  .  Both  the  Wigner 
distribution  and  the  spectrogram  are  prominent  members  of  Cohen’s  class.  The  WD,  and  its  widowed 
approximation,  is  obtained  from  (2.1)  by  setting  the  kernel  to  be  an  impulse  in  time,  t.  TFDs  have  been 
successfully  applied  to  areas  where  signals  are  localizable  in  the  time  frequency  domain  and  have  fixed 
distinct  signatures  that  permit  their  classification  and  separation.  Many  of  these  applications  are  discussed 
in  the  book  by  Cohen  [13]  and  also  in  the  book  by  Qian  and  Chen  [25].  One  successful  offering  of  TFD  is 
in  the  area  of  instantaneous  frequency  estimation  for  FM  signals  [9].  TFD  can  be  used  to  provide  an 
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accurate  estimate  of  a  jammer  uniquely  characterized  by  its  instantaneous  frequency.  The  distinctions 
between  the  FM  and  the  Gold  code  time-frequency  signatures  as  well  as  the  relative  high  strength  of  the 
jammer  and  weak  GPS  signal  power,  simplify  the  jammer  IF  estimation.  Figure  3  shows  the  Wigner 
distribution  of  a  chiip  signal  in  noise  with  SNR  =  -20  dB.  Added  to  the  chiip  signal,  a  spread  spectrum 
signal  of  JSR  =  40  dB.  The  chirp  t-f  signature  is  clear  in  the  figure,  and  its  parameters  can  be  easily 
identified  and  used  for  subspace  construction.  We  maintain  that  the  proposed  technique  for  GPS  FM 
jammer  suppression  does  not  require  the  application  of  any  particular  IF  estimator.  However,  TFD  has 
been  shown  to  outperform  other  estimators,  specifically  for  signals  with  rapidly-changing  frequency 
characteristics  [9],  with  increased  computational  cost. 

3.  Subspace  Projection  Array  Processing 

The  concept  of  subspace  projection  for  instantaneously  narrowband  jammer  suppression  is  to  remove  the 
jammer  components  by  projecting  the  received  data  onto  the  subspace  that  is  orthogonal  to  the  jammer 
subspace,  as  illustrated  in  Figure  4. 

In  GPS,  the  PN  sequence  of  length  L  (1023)  repeats  itself  Q  (20)  times  within  one  symbol  of  the  50 
bps  navigation  data.  We  use  discrete-time  form,  where  all  the  signals  are  sampled  at  the  chip-rate  of  the 
C/A  code.  We  consider  an  antenna  array  of  N  sensors,  and  the  communication  channel  is  restricted  to  flat¬ 
fading.  In  the  proposed  interference  excision  approach,  the  LNQ  sensor  output  samples  are  partitioned 
into  Q  blocks,  each  of  L  chips  and  LN  samples.  The  jammer  can  be  consecutively  removed  from  the  20 
blocks  that  constitute  one  symbol.  This  is  achieved  by  projecting  the  received  data  in  each  block  on  the 
corresponding  orthogonal  subspace  of  the  jammer.  The  jammer-free  signal  is  then  correlated  with  the 
replica  PN  sequence  on  a  symbol-by-symbol  basis.  We  first  consider  the  subspace  projection  within  each 
block.  The  array  output  vector  at  the  kth  sample  is  given  by 

D 

\(k)  =  x  s  ( k)  +  xH  ( k )  +  b(k)  =  c(k)  h  +  ^  Aiui  (k)  a(.  +  b  (k)  (2.2) 

i=i 
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where  xs,  x„,  and  b  are  the  signal,  the  jammer  and  the  white  Gaussian  noise  contributions,  respectively,  h 
is  the  signal  spatial  signature,  and  cik)  is  the  spreading  PN  sequence.  The  number  of  jammers  is  D.  All 
jammers  are  considered  as  instantaneously  narrowband  FM  signals  with  constant  amplitude 
M,(k)=exp[/^,(k)].  A,  and  a,  are  the  i- th  jammer  amplitude  and  the  spatial  signature,  respectively. 

Furthermore,  we  normalize  the  channels  and  set  h  =  /V  and  a  =  N  ,  where  •  is  the  Frobenius 

ii  iif  ii  ii/-  iiiif 

norm  of  a  vector.  The  noise  vector  b(/j  is  zero-mean,  temporally  and  spatially  white  with 

E[b(k)bT(k  +  l)]  =  0 
E[b(k)bH(k  +  l)]  =  a2S(l)lN 

where  cr  is  the  noise  power,  and  I.y  is  the  NxN  identity  matrix.  Using  L  sequential  array  vector  samples 
within  the  block,  we  obtain  the  following  L/V-by- 1  vector 


X  =  [xr(l)  xr(2)  ...  xr(Z)f 
=  Xi+XH+B 

The  vector  X„  consists  of  the  D  jammer  signals,  and  is  expressed  as 

v 

x„=Iv,. 


(2.4) 


(2.5) 


with 


V,-  =  ®  a,  =  -J=[m,(1)  «,-( 2)  ...  Uj(L)]T  ®  a, 

where  <S>  denotes  the  Kronecker  product.  Therefore, 

V  =  [V,  v2  ...  VD] 

spans  the  jammer  signal  subspace,  and  its  orthogonal  subspace  projection  matrix  is  given  by 

P  =  I^-V(V"V)-1V"  =Ii7V-lvv" 


The  projection  of  the  received  signal  vector  onto  the  orthogonal  subspace  yields 

X±  =  PX  =  PX  +  PB 


(2.6) 


(2.7) 


(2.8) 


(2.9) 
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which  excises  the  jammers  entirely.  The  signal  vector  Xs  ,  assuming  symbol  “1”  is  transmitted,  can  be 
rewritten  as 

X,  =[/>(!)  P( 2)  ...  A  q  (2.10) 


where  p(n)  is  the  C/A  code  and  the  vector  q  represents  the  spatial-temporal  signature  of  the  GPS  signal. 
The  result  of  despreading  in  the  subspace  projection  based  array  system  over  one  block  is 

y  =  q"X±=q"Pq  +  q"PB  A  yl+y2  (2.11) 


where  y\  and  y2  are  the  contributions  of  the  PN  and  the  noise  sequences  to  the  decision  variable, 
respectively.  Equation  (2.1 1)  assumes  that  the  received  satellite  signal  is  aligned  with  its  receiver  replica. 
Appendix  A  shows  that  the  projection  operation  does  not  cause  a  shift  of  the  correlation  peak,  and  as 
such,  there  is  no  bias  in  the  pseudorange.  For  simplification,  we  assume  that  the  jammers  share  the  same 
period  as  the  GPS  data  symbol.  In  GPS  systems,  due  to  the  fact  that  each  satellite  is  assigned  a  fixed 
Gold  code  [3],  and  that  the  Gold  code  is  the  same  for  every  data  symbol,  y\  is  a  deterministic  value,  rather 
than  a  random  variable  assumed  in  many  spread  spectrum  applications.  The  value  of yi  is  given  by 


where 


Define 


v1=q"Pq  =  q"(IijV-^VV//)q 


/y/A 


=  Tq~q"(V,  V  ...  V„) 


tH 


xrH 

VVA 


LN - q 

N 


l  f  D  \ 


V  i= 1  J 


(2.12) 


q"V;  =  (P® h)"  («,  ® a,)  =  (/?"u,)(h"a,.) 


H..  \,.Hr 


(2.13) 


a,  =  ■ 


h"a, 

N 


(2.14) 
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as  the  spatial  cross-correlation  coefficient  between  the  signal  and  the  zth  jammer,  and 


A  = 


P?», 

Jl 


(2.15) 


Accordingly, 


1=1 


between  the 

PN  sequence  i 

l"V/  - 

% 

(  D  \ 

ii 

(N 

(N 

if 

i-ZKI  |A|2 

V  i= l  y 

■th 


(2.16) 

(2.17) 


From  the  above  equation,  it  is  clear  thatyi  is  real,  which  is  due  to  the  Flermitian  property  of  the  projection 
matrix  P.  From  the  assumptions  in  (2.3),  y2  is  complex  Gaussian  with  zero-mean.  Accordingly, 


E[y\  =  yt  =  LN 


f  d  ,  .3 

i-SKIlA-l2 

v  «=i 


(2.18) 


Var[y]  =  Var[y2]  =  E[\y2\  ] 

=  £[q"PBB"P"q]  =  qHPE[BBH  ]Pq 


=  <r2q//PPq  =  o 

2q"Pq 

(  D 

=  <J~}\  =  ct2LN 

1/ 

M2 

^  *= 1 

/ 

(2.19) 


The  above  two  equations  are  derived  for  only  one  block  of  the  GPS  signal  symbol.  Below,  we  add 
subscript  m  to  identify  y  with  block  m  (m= 1,2 ,  Q).  By  summing  all  Q  blocks,  we  obtain  the  output  of 

the  symbol-level  despreading, 


Q 

y  =  X  y» 

m= 1 


(2.20) 


Since  ym  is  Gaussian  with  zero-mean,  then y  is  also  a  zero-mean  Gaussian  random  variable.  The  decision 
variable  yr  is  the  real  part  of  y, 

yr  =  Re[  vj  (2.21) 

The  expected  value  of  ly  is 
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E[yr]  =  E[y]=^E[ym] 


=ln 


(2.22) 


D  (  .  Q  . 

=ln  q-'Yj  K-l  Z|A„,-| 

i= 1  V  m=l 


where  aw  and  pmi  are  the  spatial  and  temporal  cross-correlation  coefficients  between  the  signal  and  the  i 
jammer  over  block  m.  Since  the  changes  in  the  spatial  signatures  of  the  signals  and  jammers  are  very 
small  compared  with  the  period  of  one  GPS  symbol  (20  ms),  a,m  can  be  simplified  to  aj  for  the  last  step  in 
(2.22).  The  variance  ofy,-  is 


1  1  ~ 

2  _  1  2  1  2 

(7  v  —  (7  v  —  /  (7 

>r  2  y  2  >'“ 


=\v2ln  e-Z  K- |2ZK,| 

^  i= 1 1  m= 1 


(2.23) 


Therefore,  the  receiver  SINR  expression  [17]  after  projection  and  despreading  is  given  by 


SINR  = 


Var[y ] 

2LN  Q-ft\a, ftKf 
i= 1  V  m=l 


(2.24) 


The  temporal  and  spatial  coefficients  appear  as  multiplicative  products  in  (2.24).  This  implies  that  the 
spatial  and  temporal  signatures  play  equivalent  roles  in  the  receiver  performance.  In  the  absence  of  the 
jammers,  no  excision  is  necessary,  and  the  SINR  of  the  receiver  output  will  become  2LNQ/<X,  which 
represents  the  upper  bound  of  the  interference  suppression  performance.  Clearly,  the  term 


2ZNZ  H2IK 

z=l  V  m= 1 


(2.25) 


in  equation  (2.24)  is  the  reduction  in  the  receiver  performance  caused  by  the  proposed  interference 
suppression  technique.  It  reflects  the  energy  of  the  signal  component  that  is  in  the  jammer  subspace.  We 
note  that  if  the  jammers  and  the  DSSS  signal  are  orthogonal,  either  in  spatial  domain  («,=0)  or  in 
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temporal  domain  (/?„„=0),  the  interference  excision  is  achieved  with  no  loss  in  performance.  In  the  general 
case,  for  FM  interference,  takes  a  small  value,  which  is  much  smaller  than  a,.  Therefore,  the 
difference  in  the  temporal  signatures  of  the  incoming  signals  allows  the  proposed  projection  technique  to 
excise  FM  jammers  effectively  with  only  insignificant  signal  loss.  The  spatial  cross-correlation 
coefficients  a ,  are  fractional  values  and,  as  such,  further  reduce  the  undesired  term  in  (2.25). 

Interference  suppression  using  arrays  is  improved  in  several  ways.  First,  the  employment  of  an 
antenna  array  can  lead  to  an  accurate  IF  estimation  of  the  jammers  [19].  Second,  in  comparison  to  the 
single  sensor  case  [  N  =  1 ,  a.  =  1  ], 


2  L 


e-zsi/u 


SIMR  =  — - - SVS - (2.26) 

<J 

multi-sensor  receivers,  at  minimum,  improve  SINR  by  the  array  gain.  This  is  true,  independent  of  the 
underlying  fading  channels  and  scattering  environment.  Finally,  spatial  selectivity,  highlighted  by  the 
role  of  a  is  used  to  discriminate  against  the  jammer  signal. 


4.  Effects  of  IF  Errors  on  the  Projection  Operation 


Errors  in  IF  may  occur  in  many  situations,  where  it  becomes  difficult  to  determine  the  IF  due  to  a  drop  in 
the  jammer  power,  presence  of  amplitude  modulations,  or  high  levels  of  cross-terms  in  the  t-f  domain. 
When  IF  estimation  errors  exist,  the  subspace  projection  operation  will  not  entirely  remove  the  jammer. 
The  un-excised  residual  jammer  at  the  projection  filter  output  is  often  significant,  specifically  for  high 
JSR.  With  no  specific  focus  on  any  particular  IF  estimator,  the  phase  errors  in  this  section  are  modeled  as 
a  zero-mean  Gaussian  white  noise  process,  motivated  by  the  fact  that  phase  errors,  directly  obtained  from 
the  analytic  signal  of  FM  in  complex  Gaussian  additive  noise,  have  wrapped  Gaussian  distributions  [20]. 
For  high  jammer  power,  the  distribution  variance  becomes  very  small  and  the  phase  errors  assume  a 
Gaussian  distribution.  Consider  a  single  jammer  with  an  estimated  unit  vector  represented  as 
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_  1  rc./Wl)+A(l))  em 2)+A(2)) 

V^L 


emL)+A(L )) 


(2.27) 


The  phase  estimation  errors  A(i)  at  different  chips  are  assumed  to  be  i.i.d  random  variables  with  a  zero 
mean  Gaussian  distribution  and  variance  a]  .  The  variance  a]  is  assumed  to  be  sufficiently  small  such  that 
most  errors  lie  inside  the  interval  [-n,  rc\.  The  projection  matrix,  constructed  from  the  inaccurate  jammer 
vector,  is 

>  =  IiiV-^UUw  (2.28) 


where  U  =  u®a  .  In  this  case,  the  output  of  the  correlator  in  one  block  is 

9  =  XfPXs+XfPW  +  XfPX„  A  yl+y2  +  y3 


(2.29) 


where  y, ,  y2  and  y3  represent,  respectively,  the  contributions  of  the  spreading  code,  the  noise  sequence, 
and  the  interfering  signal.  Due  to  phase  estimation  errors,  these  three  variables  are  random  which  renders 
equation  (2.29)  different  from  its  deterministic  counterpart  equation  (2.11).  Since  the  projection  matrix 

is  Hermitian,  yt  is  always  real.  The  mean  value  of  y{  is 


E[y  ,]  =  E[\”P\s]  =  LN  -  ±E[X?  UU"XJ 


(2.30) 


Define 


N 


(2.31) 


as  the  spatial  cross-correlation  coefficient  between  the  signal  and  the  jammer,  and 


(2.32) 


as  the  exact  and  estimated  temporal  cross-correlation  coefficient  between  the  PN  sequence  and  the 
jammer  vector,  respectively.  Therefore, 

xf  U  =  (p  <E>  h)ff(u  ®  a)  =  (p^uKh^a)  =  -JZ NajB  (2.33) 


It  can  be  readily  shown  that 
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(2.34) 


Hence, 


£[uuff]  =  e  a‘u  u7/  +  -^-(1  —  e  <J‘  )IP 

EM  =  LN(  1  -  |a|2  £[|^|2])  =  LN(  1  - 1«|2  £[P^E  ]) 

(  _  2 
,  ,2  2  ,  ,2  I  —  £  £ 

=  LN  l-\a\  <TCT*|^  + — - - 


(2.35) 


(2.36) 


(2.37) 


From  (2.35)  -  (2.37),  the  mean  value  of  y  can  be  calculated  by  the  sum 

E[y]  =  E[yl\  +  E[y2\  +  E[y3\  (2.38) 

The  mean  square  value  due  to  the  signal  is 

E[\y^]  =  E[X"PXsX?PXs] 

=  E[L2N2  -  2LXf  uu"x,  +  — ^xf  uu"x,xf  UUffXJ 
N 

=  L2N 2  -  27V2z|a|~  /?[pruuHp]  +  N2  |a|4  £[p7  uu^pp'  uuHp]  (2.39) 

=  L2N2  -  2N2L \af  \pf  + 1  -  e ^ )  + 

N2  \a\4  (2  +  4 L  |/?|2  )(1  -  )  +  L2 1/?|4 

In  deriving  the  above  expression,  p'  u  is  approximated  by  a  complex  Gaussian  random  variable,  invoking 
the  Central  Limit  Theorem.  From  (2.35)  and  (2.39),  the  variance  of  Vi  can  be  computed  as 

&»=£[*]-£2[p,]  (2-40) 

The  mean  square  values  of  y2  and  >’3  are  given  by 


53 


E[  y2  2]  =  fry2  =  £[X? PWWHPX  J 


=  E[E[X” PWWHPXS  \  P]] 
=  <72P[Xf  PPXJ 
=  a2E[X“PXs]  =  a2E[yi] 


■-  <j2LN 


1-  «' 


>1 


2  1  —  e 

H - 


P[|93|2]  =  P[Xf  PXsXf  PXJ  =  A2LE[VhPXsX?PV] 

=  A2LE[VhX  xfu  -  —  u^uu^x  xfu 

s  N 

- — uwx  xf uu"u  +  -^u"uuffx  xf UU"U] 

N  N- 

=  A2L{N2L  |a|2  \pf  -  N2  |a|2  y[ZfiE[ u^uu^p] 

-  N2  |a|2  JPfi*E[pHuuHu]  +  N2  \af  PIu^uu^pp^uu^u]} 
=  A2N2L  |a|2  {L \pf  -  2 L  [ e “CT*  +  (1  -  )/L ] 

+  [e~a‘  +  (1  -  e_°* )  /  L](L  \pf  e~a‘  +l-e~a‘)} 


Therefore,  the  variance  of  y3  is  given  by 


E[ 


y3 


]-T[r3] 


(2.41) 


(2.42) 


(2.43) 


The  expression  of  the  covariance  and  cross-correlation  between  the  three  components  of  the  decision 
variable  has  no  closed-form  expression.  Extensive  numerical  computations  show  that  the  covariance 

between  y.  and  y .  ( i  ^  j),  i ,  j  =1,2,3)  assume  small  values  relative  to  the  respective  variance  values, 

and  as  such  their  contributions  to  the  overall  variance  expression  below  can  be  ignored.  The  variance  of 
the  decision  variable  can  be  approximated  by 

&»*&» +}(&»+&»)  (2-44) 


The  above  equations  are  derived  for  only  one  block  of  the  signal  symbol.  Below,  the  subscript  m  is  added 
to  identify  y  with  block  m  (m=l,  2,  ...,  Q ),  and  should  not  be  confused  with  those  used  in  (2.29).  The 


decision  variable  yr  can  be  expressed  as 


Q  ^ 

t=Rc[ZtJ 

m—\ 


(2.45) 
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Since  the  white  Gaussian  noise  sequences  and  estimation  errors  are  independent  for  different  blocks,  the 


expected  value  and  variance  of  yr  become 

Q  ^  O 

£[v,.]  =  XRe[£[>’™]]’  Fa''[>’J  =  Zvar[>™]  (2-46) 

m= 1  m= 1 

The  above  expressions  can  now  be  used  to  generate  the  desired  receiver  SINR, 

SINR  =  -^d-  (2.47) 

Var[yr\ 

5.  Simulation  Results 

In  this  section,  we  present  simulation  results  illustrating  the  performance  properties  of  the  proposed 
projection  technique. 

Figure  5  depicts  the  receiver  SINR,  given  by  equation  (2.24),  as  a  function  of  the  input  SNR.  We 

consider  two  chirp  jammers,  each  of  a  sweeping  time  equals  to  the  GPS  C/A  symbol  duration.  The 

angle-of-arrival  (AOA)  of  the  satellite  signal  and  the  jammers  are  5°,  40°,  and  60°,  respectively.  A  simple 
two-element  array  is  considered  with  half-wavelength  spacing  (A  longer  array,  or  a  higher  dimensional 
array  could  be  used,  but  would  not  draw  a  different  conclusion).  The  satellite  PN  sequence  is  the  Gold 
code  of  satellite  SV#1,  and  the  normalized  frequency  of  the  jammers  are  from  0.01  to  0.2  and  from  0.5  to 
0.3,  respectively,  and  assumed  to  be  perfectly  estimated.  (In  all  analysis  and  simulations,  the  jammer 
waveforms  are  assumed  to  occupy  part  or  the  full  band  of  the  GPS  signal  and  it  is  properly  sampled.)  The 
SINR  of  the  single  sensor  case  is  also  plotted  for  comparison.  The  array  gain  is  evident  in  Figure  5.  In 
this  example,  |«/|=0.643,  |a;|=0.340.  Over  the  20  blocks  of  the  GPS  signal  symbol,  the  temporal  cross¬ 
correlation  is  computed  for  each  jammer.  It  is  found  that  \J3j\  is  in  the  range  [0.0049-  0.0604],  whereas 
\J3i\  is  in  the  range  [0.0067-  0.0839].  These  different  values  of  the  temporal  correlation  coefficient  are  due 
to  the  fact  that  different  blocks  of  the  GPS  symbol  capture  different  segments  of  the  chirp  jammer.  With 
the  above  values,  the  term  (2.25)  is  far  less  than  2LNQ,  which  allows  SINR  (2.24)  to  be  very  close  to  its 
upper  bound,  representing  a  jammer-free  environment. 
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Next,  we  examine  the  effect  of  phase  errors  on  the  receiver  performance.  We  consider  a  periodic 
chirp  jammer  whose  period  is  equal  to  one  block  length  of  the  GPS  C/A  symbol.  The  angle  of  arrival 
(AOA)  of  the  satellite  signal  is  5°.  A  two-element  array  is  considered  with  half-wavelength  spacing.  The 
Jammer-to-Signal-Ratio  (JSR)  is  set  to  50  dB  and  SNR  equal  to  -20  dB.  Figure  6  depicts  the  simulated 
values  of  the  receiver  S1NR  vs.  the  phase  error  variance  a\ ,  which  changes  in  the  narrow  range  [0,  0.01] 
for  all  blocks.  Small  errors  in  the  phase  estimates  will  only  allow  a  small  portion  of  the  jammer  to  escape 
the  excision  process.  Flowever,  with  high  jammer-to-signal  ratio,  this  portion  can  significantly 
compromise  the  receiver  performance.  The  above  error  range  is  selected  to  demonstrate  such  an  effect, 
and  it  is  also  typical  when  using  reliable  IF  estimators.  The  AOA  of  the  jammer  signals  are  set  to  5°,  35°, 
and  65°,  respectively.  It  is  clear  from  the  figure  that,  as  the  error  variance  increases,  the  output  SINR 
decreases.  The  SINR  of  the  single  sensor  case  is  also  plotted  for  comparison.  Unlike  the  result  of  exact  IF 
estimation,  where  antenna  arrays  bring  a  constant  3  dB  array  gain,  the  receiver  SINR  in  the  presence  of 
those  errors  is  dependent  on  the  spatial  signatures  of  the  signal  and  jammer.  For  small  spatial  cross¬ 
correlation  coefficients,  the  use  of  antenna  array  allows  the  receiver  to  be  more  robust  to  the  IF  estimation 
errors.  The  relation  between  the  receiver  SINR  and  the  jammer  AOA  is  shown  in  Figure  7.  In  this  case, 
phase  error  variance  a]  in  Figure  6  was  kept  constant  at  0.01.  It  is  important  to  observe  that  the  peak  and 
the  null  in  Figure  6  correspond,  respectively,  to  the  lowest  and  highest  values  of  the  spatial  correlation 
parameter  between  the  GPS  signal  and  jammer.  This  parameter  behavior  is  shown  in  the  same  Figure  8. 

In  the  next  simulation,  we  consider  multiple  jammers.  The  angle  of  arrival  (AOA)  of  the  satellite 
signal  is  50  degrees.  The  same  satellite  code  and  array  of  the  previous  example  are  used.  The  SNR  is  set 

2 

to  -20  dB.  Figure  9  depicts  the  simulated  values  of  the  receiver  SINR  vs.  the  phase  error  variance  , 
which  changes  in  the  range  [0,  0.01]  for  all  blocks.  JSR  is  equal  to  50  dB.  Both  single  jammer  and 
multiple  jammer  cases  are  depicted.  The  AOA  of  the  jammer  signals  are  set  to  be  35°  and  65° 
respectively.  When  two  jammers  are  present,  the  jammer  with  AOA  35°  has  a  frequency  range  [0,  /r],  and 
the  one  with  AOA  65°  [0.6 n,  0.8/r].  It  is  clear  from  the  figure  that,  as  the  error  variance  increases,  the 
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output  SINR  decreases.  The  35  degree  jammer  yields  worse  performance,  since  it  possesses  higher 
correlation  with  the  desired  signal.  With  the  two  jammers  present,  the  performance  is  close  to  the  worst 
single  jammer  performance.  The  receiver  acts  on  removing  the  jammer  that  is  least  correlated  with  the 
desired  signal  using  the  spatial  degrees  of  freedom,  not  counting  much  on  the  t-f  signature  distinction 
between  its  respective  waveform  and  the  GPS  signal.  The  other  jammer  escapes  the  spatial  filtering;  in 
which  case,  the  receiver  seeks  t-f  signature  distinctions  for  its  removal.  The  errors  in  IF  prevent  complete 
suppression. 

Conclusions 

Subspace  projection  is  a  pre-correlation  technique  that  can  effectively  reject  the  wideband  interference 
effectively  when  the  jammer  has  instantaneous  narrowband  t-f  profile.  In  this  chapter,  suppression  of 
frequency  modulation  interferers  in  GPS  using  antenna  arrays  and  subspace  projection  techniques  is 
proposed.  It  is  shown  that  both  the  spatial  and  temporal  signatures  of  the  signals  impinging  on  the  multi¬ 
sensor  receiver  assume  similar  roles  in  the  receiver  performance.  The  spatial  and  temporal  signatures' 
respective  correlation  coefficients  appear  as  square  multiplicative  products  in  the  SINR  expression. 
However,  due  to  the  Gold  code  structure  and  the  length  of  the  PN,  the  differences  in  the  temporal 
characteristics  of  the  FM  jammer  and  the  C/A  code  yield  negligible  temporal  correlation  coefficients. 
These  small  values  allow  the  receiver  to  perform  very  close  to  the  no-jamming  case,  irrespective  of  the 
satellites  and  the  jammers'  angles  of  arrival.  The  fundamental  offering  of  the  array  in  the  underlying 
interference  suppression  problem  is  through  its  gain,  which  is  determined  by  the  number  of  antennas 
employed  at  the  receiver. 

Since  the  subspace  projection  matrix  is  solely  dependent  on  the  IF  estimation,  IF  estimation  errors 
will  perturb  the  projection  matrix  and  allow  part  of  the  jammer  power  to  escape  the  projection  operation. 
The  effects  of  the  IF  estimation  errors  on  SINR  performance  of  GPS  receiver  using  array  subspace 
projection  in  the  presence  has  been  analyzed.  The  phase  errors  are  modeled  as  zero-mean  white  Gaussian, 
and  independent  over  different  chips.  The  spatial  signature  is  assumed  to  be  accurately  estimated.  The 
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analysis  and  simulation  shows  that,  although  the  IF  estimation  errors  can  affect  the  receiver  SINR 
significantly,  the  combination  of  temporal  and  spatial  signature  can  provide  more  robustness  in  the 
presence  of  IF  estimation  errors,  and  as  such,  render  better  performance  than  the  single  antenna  scheme. 


58 


Appendix  A 

For  a  single  jammer  and  a  single  antenna  receiver,  the  projected  satellite  signal,  upon  correlation  with  the 
receiver  C/A  code  of  the  same  satellite  provides 

y  =  pH  \_I  -  uuH~\p  =  pH p  ~(pH u)(uH p)  (A.l) 

Using  the  definition  of  the  cross-correlation  given  by  (2.15), 

y  =  L-Lj3j3*=L(l-\j3\2)  (A.l) 

The  above  equation  assumes  the  alignment  of  the  received  C/A  code  with  the  receiver  C/A  code.  In 
general,  this  is  not  the  case.  Define  pA  as  the  vector  containing  the  shifted  C/A  code  such  that  the  vector 
elements 


Pa  (n)  =  p(n  -A)  (A.3) 

In  this  case,  the  correlation  between  the  two  non-aligned  codes  is 

vA  =PaP-(pau)(1‘hp) 

Using  the  Gold  code  property, 

PaP  =  y  /  =  ±i,o 

and 


yA=r-W 


=  -L 


\ 

) 


For  unbiased  peak  location,  we  require 
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(A.5) 
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Since  |/?|  □  1 ,  then 
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A  tighter  bound  is 
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(A.6) 


Since  the  maximum  value  of  the  cross-correlation  is  1,  and 


L-\y 

L 


« 1  for  L 


1023,  then  (A.6)  is 


readily  satisfied. 
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Appendix  B 

With  presence  of  multiple  jammers,  the  parameters  in  (2.31)  and  (2.32)  will  assume  different  values.  For 
the  z-th  jammer,  the  spatial  cross-correlation  coefficient,  exact  and  estimated  temporal  cross-correlation 


coefficients  are  given  by 


Pi  41  ’  Pi  41 


Accordingly, 


X"u,=  (P  ®  b)"(u,®  a,)  =(pru,)(h"a, )  =^LNafi 


It  is  straight  forward  to  show  that  for  M,  (M  =  D),  jammers,  the  mean  values  of  three  terms  constituting 
the  decision  variable  are  give  as  follows: 

M  2 

£Ul]  =  iJV(l-£|a,|2£[i»,  ]) 

(B. 

M  2  1  _ 

=  LN[\  -  £ \a, |=  (e*  +  — — )] 


E[y2]  =  Q 


1V1 

£nj=ii[Xfra,]=£^Vl£[XfVU,] 

1=1 

M  1  M 

=Z4  -iLBtfV, -2JxfU„U,"U,] 

1=1  ™  m= 1 

M  1  M 

= YaJUnJLzA  IpAaM) 

1=1  ^  m-\ 

M  M  1  _  -o? 

^ALNiaJl  -afie*  -^jxma^  — - — ft) 

7=1  771=1  ^ 


where 
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is  the  spatial  cross-correlation  coefficient  between  the  mth  and  the  /"'  jammer  components.  The  mean 


square  value  of  yl  give  by 


E[\y,\  ]  =  £[xf px^xf pxj 


1 


=  vn1  -  e[2lx“  uir  x  j + — £[x;  uir  xx;  uir  xs  ] 

=  L2N~  -yl2  +  yn 


(B-6) 


y  12  =  Z2^2 H  pWIp  =Yu1LNl w  +i-e-CT' ) 


The  mean  square  values  of  the  other  two  terms  are  given  by 

&L=^[|T2|2]  =  ^[XfPWW^PXJ 

=  £[£[XfPWW"PXs  I  P]]  =  cr2is[Xf  PPXJ 


=  a-E[XyXs]  =  azE[yi\ 

M 

=  CT2IAf[l-£|a,f(e-'7?  ft, 


2  \  —  e 
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(B.7) 
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where  v32 ,  v33  and  _y34  can  be  shown  to  be: 
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Therefore,  the  variance  of  v3  is  given  by 
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It  can  be  shown  that  the  covariance  between  y.  and  y  .  ( /  ^  j  ),  i,  j  =1,2,3)  assume  small  values 

relative  to  the  respective  variance  values.  The  variance  of  the  decision  variable  can  then  be  approximated 
by 

0\vy  =  9*  +  ^  (0\v2  +  0>3 )  (B.  1 3) 

Similar  to  the  single-jammer  case,  the  overall  values  of  the  mean  and  the  variance  should  be  averaged 
over  Q  blocks. 
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Figure  3.  Wigner  distribution  for  a  periodic  PN  sequence  and  chirp  jammer  in  noise  (JSR=40dB, 

SNR=20dB). 


Figure  4.  Jammer  suppression  by  subspace  projection. 
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Figure  5:  Output  SINR  vs.  SNR. 


IF  estimation  e  rro  r  va  ria  n  c  e 


Figure  6.  Receiver  SINR  v.  error  variance 
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Figure  8.  Spatial  cross-correlation  vs.  jammer  AOA 
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Phase  estimation  error  variance 


Figure  9.  Receiver  SINR  vs.  phase  estimation  error  variance  for  multiple  jammers. 
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Chapter  3 

Array  Processing  for  Nonstationary  Interference 
Suppression  in  DS/SS  Communications  Using  Subspace 

Projection  Techniques 

1.  Introduction 

There  are  several  methods  that  have  been  proposed  for  interference  suppression  in  DS/SS 
communications,  most  have  been  related  to  one  domain  of  operation  [1],  [2],  These  methods 
include  the  narrowband  interference  wavefonn  estimation  [3],  [4],  frequency  domain 
interference  excision  [5],  zero-forcing  techniques  [6],  adaptive  subspace-based  techniques  [7], 
[8],  and  minimum-mean- square  error  (MMSE)  interference  mitigation  techniques  [9]. 

Nonstationary  interferers,  which  have  model  parameters  that  change  with  time,  are 
particularly  troublesome  due  to  the  inability  of  a  single  domain  mitigation  algorithm  to 
adequately  remove  their  effects.  The  recent  development  of  the  quadratic  time-frequency 
distributions  (TFDs)  for  improved  signal  power  localization  in  the  time-frequency  plane  has 
motivated  several  new  approaches  for  excision  of  interference  with  rapidly  time-varying 
frequency  characteristics  in  the  DS/SS  communication  systems.  Comprehensive  summary  of 
TFD-based  interference  excision  is  given  in  reference  [10],  The  two  basic  methods  for  time- 
frequency  excision  are  based  on  notch  filtering  and  subspace  projections.  Utilization  of  the 
interference  instantaneous  frequency  (IF),  as  obtained  via  TFDs,  to  design  an  open  loop  adaptive 
notch  filter  in  the  temporal  domain,  has  been  thoroughly  discussed  in  [11],  [12].  Subspace 
projection  methods,  commonly  used  for  mitigating  narrowband  interference  [13],  [14],  have  been 
recently  introduced  for  suppression  of  frequency  modulated  (FM)  interference  and  shown  to 
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properly  handle  multi-component  interference,  reduce  the  self-noise,  and  improve  the  receiver 
performance  beyond  that  offered  by  other  time-frequency  based  techniques  [15],  [16],  [17]. 

The  main  purpose  of  this  chapter  is  to  integrate  spatial  and  temporal  processing  for 
suppression  of  nonstationary  interferers  in  DS/SS  systems.  Specifically,  we  extend  the 
projection-based  interference  mitigation  techniques  in  [15],  [16],  [17]  to  multi-sensor  array 
receivers.  The  proposed  multi-sensor  interference  excision  technique  builds  on  the  offerings  of 
quadratic  time-frequency  distributions  for  estimation  of  1)  the  time-frequency  subspace  and 
time-frequency  signature  of  nonstationary  signals,  and  2)  the  spatial  signature  of  nonstationary 
sources  using  direction  finding  and  blind  source  separations.  With  the  knowledge  of  the  time- 
frequency  and  spatial  signatures,  the  objective  is  to  effectively  suppress  strong  nonstationary 
interferers  with  few  array  sensors.  The  proposed  technique  does  not  require  the  knowledge  of  the 
array  response  or  channel  estimation  of  the  DS/SS  signal,  but  it  utilizes  the  distinction  in  both  of 
its  spatial-  and  time-frequency  signatures  from  those  of  the  interferers  that  impinge  on  the  array. 
With  the  combined  spatial-time-frequency  signatures,  the  projection  of  the  data  vector  onto  the 
subspace  orthogonal  to  that  of  the  interferers  leads  to  improved  receiver  perfonnance  over  that 
obtained  using  the  subspace  projection  in  the  single-sensor  case. 

The  rest  of  the  chapter  is  organized  as  follows.  In  Section  2,  the  signal  model  is  described. 
Section  3  briefly  reviews  the  subspace  projection  technique.  We  present  in  Section  4  blind 
beamforming  based  on  subspace  projection  and  derive  the  receiver  output  signal-to-interference- 
plus-noise  ratio  (SINR).  Several  numerical  results  are  given  in  Section  5.  Section  6  concludes 
this  chapter. 
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2.  Signal  Model 


In  DS/SS  communications,  each  symbol  is  spread  into  L=T/Tc  chips,  where  T  and  T(,  are, 

respectively,  the  symbol  duration  and  chip  duration.  We  use  discrete-time  form,  where  all  signal 
arrivals  are  sampled  at  the  chip-rate  of  the  DS/SS  signal.  The  symbol-rate  source  signal  is 
denoted  as  s(n),  and  the  aperiodic  binary  spreading  sequence  of  the  nth  symbol  period  is 
represented  by  c(n,l)e±  1,  /=0,1,---,Z-1.  The  chip-rate  sequence  of  the  DS/SS  signal  can  be 
expressed  as 

d(k)  =  s{n)c{n,J)  with  k  =  nL  +  l  (3.1) 

For  notation  simplicity,  we  use  c(/)  instead  of  c(n,l )  for  the  spreading  sequence. 

We  consider  an  antenna  array  of  N  sensors.  The  propagation  delay  between  antenna 
elements  is  assumed  to  be  small  relative  to  the  inverse  of  the  transmission  bandwidth,  so  that  the 
received  signal  at  the  N  sensors  are  identical  to  within  complex  constants.  The  received  signal 
vector  of  the  DS/SS  signal  at  the  array  is  expressed  by  the  product  of  the  chip-rate  sequence  d(k) 
and  its  spatial  signature  h , 

xs(k)  =  d(k)  h  (3.2) 

The  channel  is  restricted  to  flat-fading,  and  is  assumed  fixed  over  the  symbol  length,  and  as  such 
h  in  the  above  equation  is  not  a  function  of  k. 

The  array  vector  associated  with  a  total  of  U  interference  signals  is  given  by 

(3-3) 

i= 1 


74 


where  a,  is  the  array  response  to  the  /th  interferer,  uj(k) .  Without  loss  of  generality,  we  set 

||h||‘  =  N  and  |ajj,  =  N ,  i=\, 2, —,U,  where  ||[|  is  the  Frobenius  norm  of  a  vector.  The  input 

data  vector  is  the  sum  of  three  components, 

x(k)  =  xs  ( k )  +  x„  ( k )  +  b(k)  =  d(k)h_+^  ^u^k)  +  b(£)  (3.4) 

(=1 

where  b (k)  is  the  additive  noise  vector.  In  regards  to  the  above  equation,  we  make  the  following 
assumptions. 

Al)  The  information  symbols  s{n)  is  a  wide-sense  stationary  process  with  E[s(n)s* (n)j  =  1 , 
* 

where  the  superscript  denotes  complex  conjugation.  The  spreading  sequence  c(k )  is  a  binary 
random  sequence  with  E[c(k)c(k+l)]  =  S(l),  where  S(l)  is  the  delta  function. 1 
A2)  The  noise  vector  b (k)  is  zero-mean,  temporally  and  spatially  white  with 

E  [b(k)br  (k  +  /)]  =  0,  for  all/ 
and 

E[b(k)bH(k  +  l)]  =  aS(l)IN 
T  H 

where  cr  is  the  noise  power,  the  superscripts  and  denote  transpose  and  conjugate  transpose, 
respectively,  and  I#  is  the  NxN  identity  matrix. 

A3)  The  signal  and  noise  are  statistically  uncorrelated. 


i 


This  assumption  is  most  suitable  for  military  applications  and  P-code  GPS. 
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3.  Subspace  Projection 


The  aim  of  subspace  projection  techniques  is  to  remove  the  interference  components  before 
despreading  by  projecting  the  input  data  on  the  subspace  orthogonal  to  the  interference  subspace, 
as  illustrated  in  Figure  1.  The  receiver  block  diagram  is  shown  in  Figure  2. 

A  nonstationary  interference,  such  as  an  FM  signal,  often  shares  the  same  bandwidth  with  the 
DS/SS  signal  and  noise.  As  such,  for  a  chirp  signal  or  a  signal  with  high-order  frequency  laws, 
the  signal  spectrum  may  span  the  entire  frequency  band,  and  the  sample  data  matrix  loses  its 
complex  exponential  structure  responsible  for  its  singularity.  Therefore,  the  interference 
subspace  can  no  longer  be  obtained  from  the  eigendecomposition  of  the  sample  data  matrix  [13], 
[15]  or  the  data  covariance  matrix  [14],  as  it  is  typically  the  case  in  stationary  environments.  The 
nonstationary  interference  subspace,  however,  may  be  constructed  using  the  interference  time- 
frequency  signature.  Methods  for  estimating  the  instantaneous  frequency,  instantaneous 
bandwidth,  and  more  generally,  a  time-frequency  subspace,  based  on  the  signal  time-frequency 
localization  properties  are,  respectively,  discussed  in  references  [18],  [19],  [15]. 

For  the  general  class  of  FM  signals,  and  providing  that  interference  suppression  is  performed 
separately  over  the  different  data  symbols,  the  interference  subspace  is  one-dimensional  in  an  L- 
dimensional  space.  We  note  that  since  an  FM  interference  has  a  constant  amplitude,  its 
respective  data  vector  can  be  detennined  from  the  IF  up  to  a  complex  multiplication  factor.  The 
unit  norm  nonnalization  of  this  vector  represents  the  one-dimensional  interference  subspace 
basis  vector.  Among  candidate  methods  of  IF  estimation  is  the  one  based  on  the  time-frequency 
distributions.  For  example,  the  discrete  form  of  Cohen’s  class  of  TFD  of  a  signal  x(t)  is  given  by 
[20] 
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Dtr(7,/)=^  ^  </>(m,r)x(t  +  m  +  T)x*(t  +  m-T)e  j4lTjr  (3.5) 

m=— co  t=— co 

where  (jKm,x)  is  a  time-frequency  kernel  that  could  be  signal-dependent.  The  TFD  concentrates 
the  interference  signal  power  around  the  IF  and  makes  it  visible  in  the  noise  and  PN  sequence 
background  [18],  [21].  It  has  been  shown  that,  for  linear  FM  signals,  Radon-Wigner  transform 
provides  improved  IF  estimates  over  the  TFD  [22].  Parametric  methods  using  autoregressive 
model  have  also  been  proposed  [23]. 

Other  nonstationary  interference  with  instantaneous  bandwidth  or  spread  in  the  time- 
frequency  domain  are  captured  in  a  higher-dimension  subspace.  In  this  case,  the  interference 
subspace  can  be  constructed  from  the  interference  localization  region  Q  in  the  time-frequency 
domain  (see,  for  example,  [15]).  The  subspace  of  interest  becomes  that  which  fills  out  the 
interference  time-frequency  region  Q  energetically,  but  has  little  or  no  energy  outside  Q. 

Interference-free  DS/SS  signals  are  obtained  by  projecting  the  received  data  vector  (in  the 
temporal  domain  processing,  the  vector  consists  of  data  samples  at  different  snapshots)  on  the 
subspace  orthogonal  to  the  interference  subspace. 

3.1.  Temporal  Processing 

In  the  single-sensor  receiver,  the  input  data  is  expressed  as 

u 

x(k)  =  xs(k)  +  xu(k)  +  b(k)  =  d(k)  +  '^jui(k)  +  b(k)  (3.6) 

1=1 

Using  L  sequential  chip-rate  samples  of  one  symbol  of  the  received  signals  at  time  index  k,  we 
obtain  the  following  input  vector 
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(3.7) 


\x{k )  x(k  —  1)  •••  x(k  +  L  +  l)]r 
=  [*,(*)  jc,(A:-1)  •••  xs(K-L  +  l)f 
+  [*,(*)  *H(*-1)  -  xH  (K-L  +  l)]r 
+[&(*)  b(k- 1)  •••  b(£  +  Z  +  l)f 

or  simply 

X(k)  =  Xs(k)  +  Xu(k)  +  B(k)  (3.8) 

We  drop  the  variable  k  for  simplicity,  with  the  understanding  that  processing  is  performed 
over  the  nth  symbol  that  starts  at  the  Mi  chip.  Then,  equation  (3.8)  becomes 

X  =  Xs+Xu+B  (3.9) 

Below,  we  relax  the  FM  condition  used  in  [13],  [16]  that  translates  to  a  single  dimension 
interference.  The  general  case  of  an  interference  occupying  higher  dimension  subspace  is 
considered.  We  assume  that  the  z'th  interferer  spans  M ■  dimensional  subspace,  defined  by  the 

orthonormal  basis  vectors  Vn,  Vi2,  •••,  F  M  ,  and  the  different  interference  subspaces  are 
disjoint.  Define 

r,~[ra  va  ■■■  vIMi]  (3.io) 

and  let  M  =  V'  M  as  the  number  of  total  dimensions  of  the  interferers.  With  L>M,  the  LxM 
matrix 

V  =  [V1  V2  •••  VLI],  nVj  =Q  for  i  ^  j  (3.11) 

is  full  rank  and  its  columns  span  the  combined  interference  subspace  J.  The  respective  projection 
matrix  is 

p  =  V(VHV)~1VH  (3.12) 

The  projection  matrix  associated  with  the  interference  orthogonal  subspace,  G,  is  then  given  by 
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p  =  iL-v(vHvy1vIi 


(3.13) 


When  applied  to  X,  matrix  P  projects  the  input  data  vector  onto  G,  and  results  in 

X±=PX  =  PXs+PB  (3.14) 

which  no  longer  includes  any  interference  component. 

The  single-sensor  receiver  implementing  subspace  projection  for  excision  of  a  single 
instantaneously  narrowband  FM  interferer  (i.e.,  U  =  1,  Mx  - 1  )  in  DS/SS  communications  is 


derived  in  [24].  The  receiver  SINR  is  shown  to  be 


SINR  = 


(Z-l)2  _  L- 1 

-  1-2 

1_iJ+CT<i_1)  uZ^i)+CT 


(3.15) 


For  typical  values  off,  (L-2)/(L-\)  «  1,  and  equation  (3.15)  can  be  simplified  as 


SINR  « 


N(L-l) 
cr  +  N/L 


(3.16) 


Compared  to  the  interference-free  environment,  where  the  receiver  SINR  is  LI  a,  nonstationary 
interference  (3.16)  is  achieved  by  reducing  the  processing  gain  by  1  and  increasing  the  noise 
power  by  the  self-noise  factor  of  1 IL. 

4.  Subspace  Projection  in  Multi-Sensor  Receiver 


In  this  section,  we  consider  nonstationary  interference  excision  in  multi-sensor  receivers  using 
subspace  projections.  We  note  that  if  the  subspace  projection  method  discussed  in  Section  III  is 
extended  to  an  A-elcmcnt  array  by  suppressing  the  interference  independently  in  each  sensor 
data  and  then  combining  the  results  by  maximum  ratio  combining  (see  Figure  3),  then  it  is 
straightforward  to  show  that  the  receiver  SINR  is  given  by 

SINR  »  N1LXX1  (3.i7) 

a  +  N/L 
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The  above  extension,  although  clearly  improves  over  (3.16),  does  not  utilize  the  potential 
difference  in  the  spatial  signatures  of  signal  arrivals,  and,  therefore,  is  inferior  to  the  receiver 
proposed  in  this  Section. 

4.1.  Spatio-Temporal  Signal  Subspace  Estimation 

To  construct  the  spatio-temporal  signal  subspace  of  the  interference  signals,  it  is  important  to 
estimate  both  the  time-frequency  signature  (or  subspace)  and  the  spatial  signature  of  each 
interferer.  The  IF  estimation  of  an  FM  interference  signal  based  on  time-frequency  distribution 
is  addressed  in  Section  III.  It  is  noteworthy  that  when  multiple  antennas  are  available,  the  TFD 
may  be  computed  at  each  sensor  data  separately  and  then  averaged  over  the  array.  This  method 
has  been  shown  in  [25]  to  improve  the  IF  estimation,  as  it  reduces  noise  and  cross  terms  that 
often  obscure  the  source  true  power  localization  in  the  time-frequency  domain. 

On  the  other  hand,  the  estimation  of  source  spatial  signature  can  be  achieved,  for  example,  by 
using  direction  finding  and  source  separation  techniques.  When  the  interference  signals  have 
clear  bearings,  methods  like  MUSIC  [26]  and  maximum  likelihood  (ML)  [27]  can  be  used  to 
estimate  the  steering  matrix  of  the  interference  signals.  These  methods  can  be  revised  to 
incorporate  the  TFD  of  the  signal  arrivals  for  improved  perfonnance  [28],  [29],  [30].  On  the 
other  hand,  in  fading  channels  where  the  steering  vector  loses  its  known  structure  due  to 
multipath,  blind  source  separation  methods  should  be  used  [31],  [32],  [33].  Since  the  interferers 
in  DS/SS  communications  often  have  relatively  high  power,  good  spatial  signature  estimation  is 
expected. 

More  conveniently,  the  spatial  signatures  can  be  simply  estimated  by  using  matched  filtering 
once  the  time-frequency  signatures  are  provided.  The  maximum  likelihood  estimator  for  the 
vector  a,  is  obtained  as 
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(3.18) 


a,-  =  V/V Yt(kMk)/  y^*(fc)x(fc) 

&=0  &=0  77 

where  u.(k)  is  the  estimated  wavefonn  of  the  /th  interferer.  It  is  noted  that  the  possible  phase 
ambiguity  in  the  waveform  estimation  of  ut{k)  does  not  affect  the  estimation  of  the  spatial 

signature.  For  slowly  varying  channels,  the  above  average  can  also  be  performed  over  multiple 
symbols  to  improve  the  estimation  accuracy. 

In  the  analysis  presented  herein,  we  assume  knowledge  of  the  interference  subspace  and  its 
angle-of-arrival  (AO A)  to  derive  the  receiver  SINR. 


4.2.  Proposed  Technique 


The  subspace  projection  problem  for  nonstationary  interference  suppression  in  DS/SS 
communications  is  now  considered  within  the  context  of  multi-sensor  array  using  N  array 
elements.  We  use  one  symbol  DS/SS  signal  duration  (i.e.,  L  chip-rate  temporal  snapshots),  and 
stack  L  discrete  observations  to  construct  an  NLx  1  vector  of  the  received  signal  sequence  in  the 
joint  spatio-temporal  domain.  In  this  case,  the  received  signal  vector  in  (3.4)  becomes 

[xr(k)  xr(^-l)  •••  x/^k  +  Z  +  l)] 

=  [xl(k)  xTs(k- 1)  -  xTs(K-L  +  l)]T 

(3.19) 

+  [*!(*)  *.r(*-l)  -  &K-L  + 1)] 

+  [br(k)  br(k-l)  •••  br(k  +  L  +  \)J 

or  simply 


X  =  X+X„+B 


(3.20) 


where  again  the  variable  k  is  dropped  for  simplicity. 
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In  (3.19),  the  interference  vector  in  the  single-sensor  problem,  given  by  (3.7),  is  extended  to  a 
higher  dimension.  With  the  inclusion  of  both  temporal  and  spatial  samples,  the  /nth  basis  of  the 
/'th  interference  becomes 

y  =y  0a.  (3.21) 

and 

V,=[vu  Vl2  ...  V.Mi]  (3.22) 

where  0  denotes  the  Kronecker  product.  The  columns  of  the  LNxM  matrix 

V  =  [V)  v2  •••  VG]  (3.23) 

spans  the  overall  interference  signal  subspace.  For  independent  spatial  signatures,  the  matrix 
rank  is  M.  The  orthogonal  projection  matrix  is  given  by 

P  =  Ii7V-V(V//V)-1V"  (3.24) 

The  projection  of  the  signal  vector  on  the  orthogonal  subspace  of  the  interferers’  yields 

X=PX  =  PXs+PB  (3.25) 

The  block  diagram  of  the  proposed  method  is  presented  in  Figure  4.  As  shown  in  the  next  section, 
effective  interference  suppression  can  be  achieved  solely  based  on  the  spatial  signatures  or  the 
time-frequency  signatures,  or  it  may  require  both  infonnation. 

4.3.  Performance  Analysis 

Below  we  consider  the  performance  of  the  multi-sensor  receiver  system  implementing  subspace 
projections.  Recall  that 

V"V  =  0  for  any  i,m  ^  j,n.  (3.26) 

and 

V"V  =  Mm  (3.27) 
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the  projection  matrix  P  becomes 


P  = 


-  — VV" 
N 


(3.28) 


The  signal  vector  Xs  can  be  rewritten  as 

X,=[x,r(*)  ETs(k- 1)  -  x*  +  L  +  l)]r 

=  [d(k)  hT  d(k-  l)hr  •••  d(k-L  +  \)hTJ 

=  s{n)[c{L- 1)  c(L-  2)  •••  c(0)f®h 
□  s(«)q 


where  the  LNx  1  vector 


q  =  [c(Z-l)  c(L- 2)  •••  c(0)f®hDc®h 


(3.30) 


defines  the  spatio-temporal  signature  of  the  desired  DS/SS  signal,  q  is  the  extension  of  the 
DS/SS  code  by  replicating  it  with  weights  defined  by  the  signal  spatial  signature. 

By  performing  despreading  and  beamforming,  the  symbol-rate  decision  variable  is  given  by 


yin)  =  q  Mx±(k)  =  5(«)q"Pq  +  q"PB(£)  □  yx(n)  +  y2(n) 


(3.31) 


where  yx{n)  is  the  contribution  of  the  desired  DS/SS  signal  to  the  decision  variable,  and  y2(n)  is 

the  respective  contribution  from  the  noise. 

The  SINR  of  the  array  output  becomes  (see  Appendix  A) 


f  U  2V 


SINR  = 


£:[y(Q] 
[v(*)] 


V  1=1 


J 


var 


fu  2V 


V  i=i 


cr 
H - 

N 


(  u 


(3.32) 


v  M 
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where  is  defined  in  (A. 9),  and  [);  is  the  spatial  correlation  coefficient  between  the  spatial 
signatures  h_and  a. ,  /  =  1, 2,  •  •  •  U ,  and  is  given  by 


'  N~  -1 


(3.33) 


Note  that  when  the  noise  power  is  small,  i.e.,  cr«l,  the  variance  of  becomes  dominant,  and 

the  output  SINR  reaches  the  following  upper  bound 

f  u  \2 

L-t»M 


SINR 


v  /=! 


J 


high  SNR  ,  v  \2  u 

V  1=1  J  1=1 


(3-34) 


This  result  is  affected  by  the  factors  L,  /V/;,  |p/|,  and  ,i  =  l,--,U .  On  the  other  hand,  when  the 


noise  level  is  very  high,  i.e.,  cr»l,  the  noise  variance  plays  a  key  role  in  detennining  var[v(k)], 
and  the  output  SINR  becomes 


SINR 


(  u  y 

V  i=l  J 


low  SNR  f 

a  1 


N 


\p, 


V  1=1 


N 


a 


(  u 


\P\ 

V  1=1  J 


(3.35) 


Unlike  the  high  input  SNR  case,  the  output  SINR  in  (3.36)  also  depends  on  both  N  and  cr. 
Comparing  (3.34)  and  (3.35),  it  is  clear  that  the  improvement  in  the  receiver  SINR  becomes 
more  significant  when  the  spatial  signatures  produce  small  spatial  correlation  coefficients  and 
under  high  SNR. 

Next,  we  consider  some  specific  important  cases.  When  /?.  =  0,/  =  !,•••  U,  var[ yx(nj\  =  Othe 


receiver  SINR  in  (3.32)  becomes  SINR  =  LN/a.  This  is  to  say,  the  output  SINR  is  improved  by  a 
factor  of  LN  over  the  input  signal-to-noise  ratio  (SNR)  (not  the  input  SINR!).  This  implies  that 
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the  interferers  are  suppressed  by  spatial  selectivity  of  the  array  and  their  suppression  does  not 
cause  any  distortion  of  the  temporal  characteristics  of  the  DS/SS  signal.  The  DS/SS  signal  in 
this  case  enjoys  the  array  gain  that  contributes  the  factor  N  to  the  SINR. 

For  a  single  FM  interferer  (U=  1,  My=  1),  equation  (3.32)  becomes 


SINR  « 


(3.37) 


It  is  easy  to  show  that  SINR  in  (3.37)  monotonously  decreases  as  |pj|  increases,  and  the  lower 
bound  of  the  SINR  is  reached  for  Pj=l,  which  is  the  case  of  the  desired  DS/SS  signal  and  the 
interference  signal  arriving  from  the  same  direction.  With  a  unit  value  of  pj, 


SINR  « 


N(L-l) 

N 

- her 

L 


(3.38) 


This  result  is  the  same  as  that  of  the  single-sensor  case  developed  in  [16],  except  for  the 
appearance  of  the  array  gain,  N,  for  the  desired  DS/SS  signal  over  the  noise.  This  equation  also 
coincides  with  (3.17).  That  is,  the  independent  multi-sensor  subspace  projection,  illustrated  in 
Figure  3,  results  in  the  same  output  SINR  with  the  proposed  multi-sensor  subspace  projection 
method  when  |Pj|=0. 


On  the  other  hand,  the  maximum  value  in  (3.36)  corresponds  to  J3=  0,  and  is  equal  to  SINR  = 
LN/a,  as  discussed  above.  For  the  illustration  of  the  SINR  behavior,  we  plot  in  Figure  5  the 
SINR  in  (3.36)  versus  |pj|  for  a  two-sensor  array,  where  L=64,  and  one  FM  jammer  is  considered 


with  M  =  7.  The  input  SNR  is  OdB. 
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4.4.  Remarks 


Given  the  temporal  and  spatial  signatures,  the  proposed  technique  simplifies  to  two  consecutive 
tasks.  The  first  is  to  estimate  the  spatio-temporal  signature.  When  using  multiple  antenna 
receivers,  a  basis  vector  of  the  orthogonal  projection  matrix  is  obtained  by  the  Kronecker  product 
of  a  jammer’s  temporal  signature  and  its  spatial  signature  that  results  in  the  LNxLN  orthogonal 
project  matrix  instead  of  LxL  in  the  single  antenna  case.  The  second  task  is  jammer  suppression 
via  subspace  projection.  This  involves  the  multiplication  of  an  LNxLN  matrix  and  an  LNxl 
vector. 

Note  such  increase  in  computations  is  natural  due  to  increase  of  dimensionality.  It  is 
noteworthy  that  array  processing  expands  overall  space  dimensionality  but  maintains  the  jammer 
subspace  dimension.  As  a  result,  it  yields  improved  SINR  performance  over  temporal  processing 
or  spatial  processing  only  methods. 

5.  Numerical  Results 

A  two-element  array  is  considered  with  half-wavelength  spacing.  The  DS/SS  signal  uses  random 
spreading  sequence  with  L= 64.  The  AOA  of  the  DS/SS  signal  is  0  degree  from  broadside 
(^=0°). 

We  consider  two  interference  signals.  Each  interference  signal  is  assumed  to  be  made  up  of 
uncorrelated  FM  component  with  M.  =  l,i  =  1,2.  The  overall  interference  subspace  is  M  =  14. 

The  AO  As  of  the  two  interferers  are  9/=[40°,60o].  The  respective  spatial  correlations  in  this 

example  are  10^=0.53  and  |(32|=0.21.  Note  that,  in  the  subspace  projection  method,  the  output 

SINR  is  independent  of  the  input  jammer- to-signal  ratio  (JSR),  since  the  interferers  are  entirely 
suppressed,  regardless  of  their  power.  Figure  6  shows  the  receiver  SINR  versus  the  input  SNR. 
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The  upper  bounds  correspond  to  interference-free  data.  For  high  input  SNR,  the  receiver  SINR 
is  decided  by  the  induced  signal  distortion,  described  by  the  variance  given  in  (A.  10).  It  is 
evident  from  Figure  6  that  the  two-antenna  receiver  outperforms  the  single-antenna  receiver  case 
by  a  factor  much  larger  than  the  array  gain.  Since  the  output  SINR  in  the  two-antenna  receiver 
highly  depends  on  the  spatial  correlation  coefficients,  the  curves  corresponding  to  a  two-sensor 
array  in  Figure  6  will  assume  different  values  upon  changing  Pj  ,  or/and  p9  .  The  best 

performance  is  achieved  at  P1=P2=0. 

Figure  7  shows  the  receiver  SINR  versus  the  number  of  chips  per  symbol  ( L ).  We  let  L  vary 
from  8  to  4096,  whereas  the  input  SNR  is  fixed  at  0  dB.  The  two  interference  signals  are  incident 

on  the  array  with  angles  9/=[40o,60°].  They  are  assumed  to  maintain  their  time-frequency  spread 

with  increased  value  of  L.  As  such,  the  respective  dimensions  of  their  subspaces  grow 
proportional  to  the  number  of  chips  per  symbol.  In  this  example,  the  dimension  of  each 
interference  signal  is  assumed  to  be  10  percent  of  L  (round  to  the  nearest  integer).  The  ouput 
SINR  improvement  by  performing  array  processing  at  different  L  is  evident  from  this  figure.  It 
is  seen  that,  unlike  the  case  of  the  instantaneously  narrowband  FM  interference,  where  the  output 
SINR  increases  rapidly  as  L  increases,  the  output  SINR  in  the  underlying  scenario  ceases  to 
increase  as  L  assumes  large  values.  This  is  because  the  rank  of  the  interference  signal  subspace 
increases  with  L. 

In  Figure  8  we  investigate  the  receiver  SINR  performance  versus  the  number  of  array  sensors. 
In  this  figure,  L  is  set  at  64,  and  the  input  SNR  is  0  dB.  Two  interference  signals  composed  of 
uncorrelated  FM  components  are  considered,  and  M.  =7,/  =  1,2,  are  assumed.  Two  examples 

are  used  to  examine  the  effect  of  different  AOAs.  In  the  first  example,  0y=[4O°,6O°].  The  output 
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SINR  improves  sharply  as  the  number  of  array  sensors  increases  from  one  to  three,  beyond 
which  the  improvement  becomes  insignificant.  The  differences  in  the  above  AOAs  of  the 
desired  DS/SS  signal  and  the  interference  signals  are  relatively  large,  and  a  small  number  of 
array  sensors  leads  to  negligible  spatial  correlation  coefficients.  We  also  show  a  case  with 

closely  spaced  interference  signals  where  0y=[5°,15°].  In  this  case,  the  output  SINR  slowly 

improves  as  the  number  of  array  sensors  increases. 

It  is  noted  that,  when  we  consider  a  specific  case,  the  output  SINR  does  not  increase 
monotonously  with  the  number  of  array  sensors.  This  is  because  the  relationship  between  the 
spatial  correlation  coefficient  and  the  AOAs  is  by  itself  not  monotonous.  Nevertheless,  when  we 
consider  the  general  case  with  different  AOA  combinations,  high  number  of  array  sensors  often 
reduce  the  spatial  correlation  coefficients. 

6.  Conclusions 

In  this  chapter,  subspace  projection  techniques  were  employed  to  suppress  nonstationary 
interferers  in  direct  sequence  spread  spectrum  (DS/SS)  communication  systems.  Interference 
suppression  is  based  on  the  knowledge  of  both  the  interference  time-frequency  and  spatial 
signatures.  While  the  former  is  based  on  instantaneous  frequency  information  that  can  be  gained 
using  several  methods,  including  time-frequency  distributions,  the  later  can  be  provided  from 
applying  higher  resolution  methods  or  blind  source  separation  techniques  to  the  signal  arrivals. 

The  differences  between  the  DS/SS  signal  and  interference  signatures  both  in  the  time- 
frequency  and  spatial  domains  equip  the  projection  techniques  with  the  ability  to  remove  the 
interference  with  a  minimum  distortion  of  the  desired  signal. 
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The  receiver  performance  based  on  subspace  projections  was  analyzed.  It  was  shown  that  the 
lower  performance  bound  is  obtained  when  the  sources  have  the  same  angular  position.  In  this 
case,  the  problem  becomes  equivalent  to  a  single-antenna  receiver  with  only  the  presence  of  the 
array  gain.  On  the  other  hand,  the  upper  bound  on  performance  is  reached  in  the  interference- 
free  environment  and  also  corresponds  to  the  case  in  which  the  spatial  signature  of  the 
interference  is  orthogonal  to  that  of  the  DS/SS  signal. 

Numerical  results  were  presented  to  illustrate  the  receiver  SINR  dependency  on  spatial 
correlation  coefficient,  input  SNR,  and  the  PN  sequence  length. 
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Appendix  A 

To  derive  the  output  SINR  expression,  we  use  s(n)  =  +1  (the  output  SINR  is  independent  of  s(n) 
and  same  result  follows  when  s(n)  =  -1).  Then, 


E[yi(n)\  =  E[  qffPq] 


=  E  q"  I--W"  q 
l  N  J 


£[q"q]-T£[q"vv"q] 


u  =1  M,  =1  U  =1  M-)  =1 


It  is  straightforward  to  show  that 


q*VM=(c®h)"  (^,®a,)  =  (crr,„)®(h“a,)  =  »JS/£^,(0c(0 


Using  the  orthogonal  property  of  the  spreading  sequence  Al),  (A.l)  becomes 


U  mh  L- 1 


U  L- 1 


£[y,(n)]  =LN-NE  £  A  X  X WMOX#  X  X  K*  ft )cft) 

z'j=l  Wi=l/i=0  /2=1  m2  =14=0 


U  M<1  L- 1 

I  „  |2  I 


=lat-jv2I/?|  XXMO  AO 


=1  m= 1  /=0 


Due  to  the  zero-mean  property  of  noise  (assumption  A2)  E[  y2(n)]  =  0 .  Accordingly, 


E[y(n)\  =  E[yi(n)\  =  N 


It  is  clear  from  (A. 4)  that  the  increase  in  the  space  dimensionality  from  L  to  NL  does  not  simply 
translate  into  a  corresponding  increase  in  the  desired  mean  value,  or  subsequently  in  the 
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processing  gain.  Also,  from  assumption  A3),  the  cross-correlation  between  yt  (n)  and  y2(n)  is 


E [  v* (n)y2 («)]  =E[yx  ( n)y2 (»)]  =  0  ( A.5) 

Therefore,  the  mean  square  value  of  the  decision  variable  is  made  up  of  only  two  terms, 


\y(nf  =E  \yx(n)\2  +E  \ y2(n) 


(A.  6) 


The  first  term  is  the  mean  square  value  of  yt  (n) .  From  (3.26),  we  have 
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where 


and 


Mj  Mx  L- 1 

^  =  Zr>,w,;(o 

m,  =1  m-)  =1 L  =0 


(A.  8) 
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(A.9) 


In  practice,  y,  takes  negligible  values,  and  equation  (A.7)  can  be  simplified  to 
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The  value  of  ^  depends  on  the  type  of  interference  signals.  Specifically,  when  the  it h 


interference  signal  is  made  up  of  a  single  FM  or  a  number  of  uncorrelated  FM  signal  components, 
then  the  basis  vectors  are  of  constant  modulus,  and 


6  = 


L 


(A.  11) 


The  second  term  of  (A. 6)  is  the  mean-square  value  of  y2(n) , 


|j;2(n)|2]  =  JE[q"PB(A:)B^(A:)Pffq] 

=  ctE [q^PP^q]  =  aE[ q"Pq]  =  ctn[ L-fjMi  \ \ 

V  i=l 


The  variance  of  y(n)  is  given  by 

var [y{n)\  =  e\ \y{nj[~\  -E2  [yin)] 
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(A.  13) 


Equation  (3.32)  follows  by  using  the  results  of  (A. 4)  and  (A.  13). 
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Figure  1.  Jammer  suppression  by  subspace  projection. 
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Figure  2.  Block  diagram  of  single-sensor  subspace  projection. 
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Figure  3.  Block  diagram  of  independent  multi-sensor  subspace  projection. 
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Figure  4.  Block  diagram  of  proposed  multi-sensor  subspace  projection. 
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Figure  5.  Output  SINR  versus  |Pj|  (input  SNR=0dB,  L= 64,  U=  1,  Af=7). 
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Figure  6.  Output  SINR  versus  input  SNR  (L=64,  U=2,  M1=M2=7,  qD=0°,  qj=[40°,60°]). 
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Figure  7.  Output  SINR  versus  the  number  of  chips  per  symbol  ( L )  (input  SNR=0dB,  U= 2, 

Mj=M2=7,  00=0°,  0 ^=[40°, 60°]). 
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Figure  8.  Output  SINR  versus  the  number  of  array  sensors  (input  SNR=0dB,  Z,=64,  U=  2,  My=M\=l 


Chapter  4 

Performance  Analysis  of  GPS  Receivers  in  Impulsive  Noise 

1.  Introduction 

Impulsive  noise  is  encountered  in  many  environments.  It  has  been  shown  that  this  type  of  noise  has  a 
significant  effect  on  satellite-mobile  radio  systems  [1],  One  major  source  of  impulsive  noise  is  the 
automotive  ignition  systems.  The  frequency  range  of  automotive  ignition  noise  is  100  MHz  ~  10  GHz  [2], 
which  extends  over  the  GPS  carrier  frequency  of  1.58  GHz  and  its  1  MHz  bandwidth.  The  impulsive 
noise  from  the  ignition  systems  shows  a  random  nature  of  the  amplitudes,  inter-arrival  times,  and 
durations  of  noise  bursts.  This  noise  might  be  a  hidden  menace  to  GPS  receivers  since  the  ignition  wire 
and  the  body  of  the  automobile  act  as  a  radiation  antenna.  Another  source  of  impulsive  signals  is 
ultrawideband  (UWB)  signals  that  cover  the  GPS  operating  band  and  find  increasing  number  of 
applications  in  outdoor  and  indoor  environments. 

This  chapter  considers  the  effect  of  impulsive  noise  on  the  GPS  receiver  performance,  specifically,  its 
delay  lock  loop  (DLL).  We  use  Middleton  noise  [3]  and  generalized  Cauchy  noise  [4],  which  are  two  of 
the  canonical  models  that  are  widely  used  to  describe  the  statistics  of  the  impulsive  noise.  We  also 
examine  the  impact  on  the  GPS  receiver  by  automotive  ignition  noise  and  UWB  noise  models  based  on 
experiment  data  provided  in  [5],  [6],  The  central  limit  theorem  (CLT)  is  applied  to  characterize  the  noise 
components  resulting  from  the  correlation  with  the  early,  late,  and  punctual  C/A  reference  code.  It  is 
shown  that  the  correlator  noise  is  also  Gaussian.  The  variance  of  the  discriminator  is  derived  and  shown 
to  be  highly  dependent  on  the  signal-to-noise  ratio  (SNR),  sample  rate,  and  precorrelation  bandwidth. 
Computer  simulations  are  performed  and  compared  with  theoretical  results. 
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2.  The  Impulsive  Noise  Models 


2.1.  Middleton  noise  model 


Middleton  impulsive  noise  model  [3]  is  composed  of  both  Gaussian  and  impulsive  noise  components. 
The  probability  density  function  (PDF)  is  defined  by  two  parameters  A  and  r  .  The  term  A  represents  the 
product  of  the  average  duration  of  the  pulse  and  the  average  number  of  pulses  occurring  in  unit  time. 
Small  values  of  A  increases  the  impulsiveness  of  the  noise,  whereas  large  values  of  A  move  the  model 
closer  to  a  Gaussian  distribution.  The  parameter,  r  =oG2  / a,2 ,  is  defined  as  the  Gaussian-to-impulsive 


power  ratio,  where  oG2  represents  the  Gaussian  noise  power  and  o,2  is  the  impulsive  noise  power.  The 
total  noise  power  is 

o2=oG2+o,2  (4.1) 


The  PDF  is  given  by 


P(z)  =  % 
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where 


i  =  (m,A)  +  r 

m  1  +  r 

If  A  is  sufficiently  small,  then  we  can  simplify  the  model  by  only  keeping  the  dominant  terms  in  (4.2) 
corresponding  to  m  =  0,  1,  2  [7],  In  this  case,  the  noise  PDF  function,  shown  in  Figure  1(a),  can  be 
approximated  by 
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(4.4) 


A  sample  of  Middleton  noise  sequence  is  shown  in  Figure  1(b). 
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2.2.  Generalized  Cauchy  noise  model 


A  generalized  Cauchy  probability  density  function  is  defined  in  terms  of  three  parameters,  a2 ,  k  >  0  and 
»>0[4], 


P(z)  = 


B 

{i+i[M]k}-i/k 

u  A 


A  =  ra2  r(1/k)l1/2  B_kn-'/kr(n  +  l/k) 
T(3/k)  ’  2AT(u)T(l/k) 


In  the  above  equation,  T(-)  is  the  Gamma  function,  given  by 


(4.5) 


T(a)  =  J  xale~xdx  (4.6) 

0 

The  parameter  k  controls  the  impulsiveness  of  the  noise,  whereas  v  controls  the  noise  variance  and  a2  is  a 
scale  parameter.  In  the  example  shown  in  Figure  2,  we  set  k  =  0.2,  v  =  40,  and  a2  =  4.7. 


2.3  Ignition  noise  model 

The  motor  ignition  noise  based  on  the  experiment  data  [5]  is  a  statistical  process  consisting  of  two 
distributions:  a  Weibull  distribution,  due  to  the  high  power  peaks,  and  a  Gaussian  distribution,  due  to 
other  low  power  values.  These  distributions  are  thown  in  Figures  3(a)  and  3(b),  respectively.  The  Weibull 
distribution  PDF  is  given  by 

p(z)  =  atTaza-V(z/b)2  (4.7) 

where  a  =  1.14  and  b  =  4.00.  The  inter-arrival  times  between  successive  noise  peaks  are  shown  to  be 
exponentially  distributed,  as  depicted  in  Figure  3(c).  The  sample  noise  sequence  is  provided  in  Figure 
3(d). 

2.4.  UWB  noise  model 

The  UWB  type  of  noise  generated  using  multiple  UWB  signal  sources  is  reported  in  [6]  based  on  several 
data  measurements.  It  was  concluded  that  the  amplitude  of  the  aggregate  UWB  signals  approximates  a 
Rayleigh  distribution,  which  is  shown  in  Figure  4(a).  Figure  4(b)  gives  the  sample  sequence  of  the  noise 
with  unit  variance. 


102 


3.  DLL  Performance  Under  Impulsive  Noise 


The  major  operations  for  code  synchronization  in  the  DLL  are  the  cross-correlations  performed  of  the 
incoming  data  and  the  receiver  reference  code.  Figure  5  shows  three  pairs  of  correlators  required  to 
produce  three  in-phase  components,  IE  ,  Ip  ,  IL  and  three  quadraphase  components  QE  ,  Qp  ,  QL  , 
respectively  corresponding  to  the  early  (E),  punctual  (P)  and  late  (L)  reference  C/A  codes.  With  the  above 
six  components,  the  receiver  could  construct  at  least  three  different  DLL  discriminators  8,  namely 
Coherent 

D  =  (IE -IL)sign(Ip)  (4.8) 


Early-minus-late  power  (noncoherent) 


D  =  (Ie2  +Qe2)-(1l2  +Ql2) 


(4-9) 


Dot-product  (noncoherent) 


D  =  (IE  -  IL)IP  +  (Qe  '  Ql  )Qp 


(4.10) 


where  sign(Ip)  is  the  sign  of  the  navigation  message  data  bit.  Ideal  synchronization  is  reached  by  finding 
the  location  of  the  correlation  peak.  Commonly,  the  discriminator  determines  the  peak  correlation 
location  by  reaching  zero  output  value.  The  output  of  correlator  j,  corresponding  to  the 
Figure  5:  GPS  DLL  cross-correlation  process  early,  punctual  or  late  stage,  when  the  summation  is 
performed  over  T  (usually  0.001)  seconds,  can  be  written  as  [9]: 

=ij  +rhj  =  V2SMR(Tj)c°s(p  +r|ij,  Qj  =  +  4^  =  V^SMRfr^sintp  +  (4.11) 

where  S  is  the  signal  power,  M  is  the  number  of  samples  used  in  cross-correlation  computations  (usually 
it  is  an  integer  multiple  of  1023),  cp  is  the  residual  phase  tracking  error  at  the  time,  R(tj)  is  the  cross¬ 
correlation  function  between  the  incoming  C/A  code  and  the  reference  code  corresponding  to  stage  j  for  a 
delay  ij5  and  and  i]0|  are  the  in-phase  and  quadraphase  noise  components  of  the  correlator  outputs. 


Clearly,  if  gi ;  represents  the  reference  code  samples,  and  n  ^  |  and  n  ()|  i  represent  the  input  noise  samples, 
then  the  noise  components  can  be  expressed  as 
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(4.12) 


n  Xa  "  >  ^=£1  ='V 


3.1.  Sample  rate 

The  C/A  code  consists  of  1023  chips.  In  order  to  generate  the  proper  early  and  late  correlation  functions 
for  the  discriminator,  the  sample  rate,  i.e.  the  number  of  samples  per  chip,  becomes  dependent  on  the 
early  and  late  correlator  spacing  d.  For  d  =  m/n  chips,  where  m  and  n  are  integers,  forming  an  irreducible 
fraction,  then  the  number  of  samples  per  chip  should  be  no  less  than  2n,  if  m  is  not  a  multiple  of  2,  or  no 
less  than  n,  if  m  is  a  multiple  of  2.  In  Figure  6,  we  choose  the  coherent  discriminator  and  set  the  early-late 
correlator  spacing  to  1  chip.  Accordingly,  the  number  of  samples  per  chip  should  be  2.  Without  multipath 
and  noise,  the  discriminator  output,  which  is  the  difference  between  the  early  and  late  correlators, 
assumes  zero  value  at  synchronization. 

3.2.  Precorrelation  filtering 

The  precorrelation  filter  in  GPS  receiver  is  used  to  suppress  out-of-band  noise  and  interference.  However, 
filtering  may  cause  correlation  changes  and  removes  the  sidelobes  of  C/A  code  spectrum.  Figure  7  shows 
the  GPS  frequency  spectrum  transformation  due  to  the  application  of  Butterworth  filter. 

The  effect  of  the  filter  on  the  C/A  code  autocorrelation  function  is  broadening  its  peak  [10],  as  shown 
in  Figure  8.  Therefore,  the  early-late  correlator  spacing  has  to  be  selected  sufficiently  large  to  keep  the 
discriminator  properly  functioning  over  a  linear  range.  As  evident  from  Figure  8,  the  correlation  function 
is  not  absolutely  symmetric,  which  leads  to  slight  deflection  of  zero-delay  error  point.  The  ideal  band 
limited  filtered  noise  has  a  sine  autocorrelation  function,  rather  than  a  delta  function.  In  baseband, 

sin(jiBfx) 


RfnOO  =  ' 


7lBfT 


(4.13) 


where  Bf  is  the  precorrelation  filter  bandwidth  of  2  MHz.  Noise  samples  spaced  by  0.5  microseconds  are 
uncorrelated.  Oversampling  will  generate  correlated  noise  samples  which  will  adversely  affect  the 
discriminator  performance.  Correlated  and  uncorrelated  noise  samples  distort  the  early  and  late 
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correlation  functions,  causing  the  discriminator  output  to  assume  non-zero  values  at  the  correct  signal 
propagation  delay.  Figure  9  shows  the  distorted  correlation  and  discriminator  function  in  -10  dB 
impulsive  noise. 

To  evaluate  the  discriminator  error,  we  use  the  dot-product  discriminator  as  an  example.  From  (4.10) 
and  (4.11),  the  discriminator  output  is: 

D  =  (Ie  - 1  l  )I  p  +(Qe -Ql)Qp  =0e  ->l  +tIie  —  hiL  X*p  +  t1ip)  +  (4e  “  9 i,  +11qe  —  4ol  X4p  +  4qp  )  (4.14) 

Thus,  the  discriminator  variance  is 

Var(D)  —  E[ip  (r|IE  —  r|IL)  +r|IP  (r|IE— r|IL)  +  qP  (r|gE  —  r|gL)  +  vj QP  (PqE  —  PqL)  ]  (4.15) 

where 


E[0liE  4il)  1  —  E[(t)qE  r|pL)  ]  —  2(E[t|ie  ]  E[r|IEr|1L]) 

2PJdM/fiM-l)(Ud-./k)Sm('f;'BT2|:(M-.)(l-./t)S'"('f;)B-)  (4.16) 


i7iBf  /Bs 


i7iBf  /Bs 


-tx'(M-i)(i-|kd-i|/k)siii(i;B:)B')] 


i7iBf  /Bs 


E[r|n>  (hiE  4il)-]  —  E[r|0P‘(riQE  u|qL )_ ] 

=  2P/[M  +  2V  (M  -  i)(l  -  l,k)5WM][dM  +  ‘y'(M  -  1X1  +  d  -i/k) 

ijiBf/Bs  i7tBf/Bs 


+  2  g  (M  -  i)(l  -  i  /  k) sin(i,lB' ' B-  >  -  (M  -  i)(l  -  |kd  -  i| ;  k)  sin(if'  )B'  >] 


brBf  /Bs 


inB f  /  Bs 


(4.17) 


The  derivation  of  the  above  equations  is  lengthy  and  given  in  Appendix  A. 

4.  Simulations 


In  this  section,  we  present  the  results  obtained  from  20,000  Monte  Carlo  trials,  generated  based  on  the 
noise  models  described  before.  The  purpose  is  to  evaluate  the  effects  of  the  precorrelation  bandwidth  and 
sample  rate  on  the  early-late  discriminator  error  variance.  We  compare  the  analytical  results  with  the 
simulations  results.  The  first  set  of  simulation  provides  the  details  of  how  the  discriminator  error  variance 
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changes  according  to  the  precorrelation  bandwidth.  The  second  set  of  simulations  shows  the  discriminator 
statistics  in  terms  of  the  sample  rate. 

In  the  first  set  of  simulations,  we  compute  the  values  of  discriminator  error  variance  over  different 
precorrelation  bandwidths  of  2,  4,  6,  8,  and  10  MHz.  The  sample  rate  is  fixed  to  30  samples  per  chip,  and 
the  early-late  correlator  spacing  is  set  to  0.4  chips.  Both  UWB  noise  and  Middleton’s  impulsive  noise  are 
considered.  The  signal-to-noise  ratio  (SNR)  is  set  to  -30  dB.  Figure  10  shows  that  the  discriminator  error 
variance  decreases  as  the  precorrelation  bandwidth  increases.  It  is  evident  that  the  analytical  and  the 
simulation  results  are  similar.  The  result  demonstrates  that  a  narrow  precorrelation  bandwidth 
compromises  the  discriminator  performance.  For  a  clear  view,  we  divide  the  discriminator  error  variance 
based  on  (4.15)  into  four  components,  (a)  E[i p2  (r) IE  —  rj IL ) 2  ]  ,  (b)  E[q P 2  ( r) QE  —  ii QL ) 2  ]  ,  (c) 

E[ii „> 2 (r) IE  —  ti il ) 2 1  ,  and  (d)  E[r|QP2(r|QE -r|QL)2]  .  Figure  11  shows  the  simulation  and  the  analytical 

behaviors  of  those  components,  which  are  clearly  similar. 

The  second  set  of  simulations  also  applies  UWB  noise  and  Middleton’s  noise.  Figure  12  shows  the 
different  values  of  discriminator  error  variance  with  different  sample  rates  of  5,  10,  15,  20,  and  30 
samples  per  chip,  under  SNR  of  -30  dB.  The  precorrelation  bandwidth  is  fixed  as  2  MHz  and  the  early- 
late  correlator  spacing  is  set  to  0.4  chips.  We  find  that  the  discriminator  error  variance  increases  as  the 
sample  rate  increases.  This  simulation  results  agree  with  the  corresponding  analytical  results.  Figure  13 
depicts  the  close  simulation  and  analytical  values  of  the  four  components,  which  sum  up  to  the 
discriminator  error  variance  of  Figure  12. 

By  examining  the  contributions  of  the  four  components  (a),  (b),  (c)  and  (d)  to  the  discriminator  error, 
we  find  that  both  (a)  and  (b)  have  the  signal  power  terms  tp"  and  qp  ,  whereas  (c)  and  (d)  only  include 
the  noise  power  terms.  This  means  that  (a)  and  (b)  dominate  the  discriminator  error  variance  at  high 
SNRs,  whereas  (c)  and  (d)  play  a  more  important  role  at  low  SNRs.  Figure  14  shows  the  values  of  the 
aforementioned  components  at  different  SNR  of  -90,  -70,  -50,  -30  and  -10  dB,  with  the  sample  rate  of  10 
samples  per  chip  and  an  early-late  correlator  spacing  of  0.4  chips. 
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5.  Conclusions 


Several  impulsive  noise  models  were  adopted  to  evaluate  their  effect  on  the  delay  lock  loop  of  the  GPS 
receiver.  We  defined  the  required  sample  rate  for  the  early-late  discriminator  and  performed  statistical 
analyses  of  discriminator  error  variance  induced  by  the  impulsive  noise,  the  precorrelation  bandwidth, 
and  the  sample  rate.  The  required  sample  rate  is  dependent  on  the  early-late  correlator  spacing.  A 
narrower  spacing  requires  more  samples  per  chip.  The  use  of  the  precorrelation  filter  may  change  the 
correlation  between  noise  samples,  which  subsequently  affects  the  discriminator  error.  Narrower 
precorrelation  bandwidth  and  higher  sample  rate  lead  to  increased  discriminator  errors.  With  reasonable 
and  careful  selections  of  the  sample  rate  and  the  precorrelation  bandwidth,  the  GPS  synchronization  error 
can  be  limited. 
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Appendix  A 


We  begin  with  the  variance  expression  (4. 15) 

E[ (t|iE  -niL)2]  =  E[nIE2  ]  +  E[ rj  1L2 ]  -  2E[r|IEr|IL] 


We  note  that, 


M 

E[  r)1E2  ]  =  E[  (^n(i)g(i))2  ] 

i=l 


(A.l) 


M  MM 

=E[  Z(n(i)g(i))2  +  2n(l)g(l)£n(i)g(i)  +  2n(2)g(2)£n(i)g(i)  +  -  +  2n(M  -  l)g(M  -  l)n(M)g(M)  ]  (A.2) 

i=l  i=2  i=3 

where  n(i)  represents  the  noise  sample,  g(i)  is  the  C/A  code  sample,  and  M  is  the  total  number  of  samples. 
For  simplification,  we  use 

Rf(i)  =  1  —  i/k  for  i  =  1,  2,  ...,  k  -  1;  Rf(i)  =  0  for  i  =  k,  k+  1,  ...,  M  (A.3) 


where  k  is  the  number  of  samples  per  chip,  to  approximate  the  correlation  function  between  the  filtered 
and  the  reference  C/A  code,  and 


Rf„(i)  =  P„ 


sin(inBf  /  Bs) 
ijiBf  /  Bs 


(A.4) 


for  the  noise  correlation  function,  in  which  Pn  is  the  noise  variance.  Thus,  ( A2)  becomes 


E[nIE2] 


k-1 

[M  +  2£  (M  —  i)(l  —  i  /  k) 


i=l 


sin(i7iBf /Bs) 
inBf/Bs 


(A.5) 


where  the  total  number  of  samples  M  =  1 023k.  Similarly, 


E[  B1L2  ]  =  E[  riIE2  ]  =  E[  r)lp2  ] 


(A.6) 


Moreover, 

M  M 

E[t|iEt|il]  =E[(£n(i)g(i))  (£n(i)g(i +  kd))  ] 

i=l  i=l 

M 

=  E[  £  n(i)g(i)n(i)g(i  +  kd)  ] 

i=l 

M  M 

+  E[ n(l)g(l)£n(i)g(i  +  kd)  +  n(2)g(2)£n(i)g(i  +  kd)  +•••  +  n(M  - l)g(M  - l)n(M)g(M  +  kd)  ] 

i=2  i=3 
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E[ n(2)g(2)n(l)g(l  +  kd)  +  n(3)g(3)  £n(i)g(i  +  kd)  +  •  ••  +  n(M)g(M)  jn(i)g(i  +  kd)  ] 


Pn{(l-d)M+  £  (M-i)(l-d-i/k)  — 


sin(kBf /B.)  ,  k+^-‘  ..  ,,  ,  sin(ijtBf  /B.) 


i7iBf  /  Bs 


-+  £  (M-i)(l-|kd-i|/k)  — 


i7iBf  /  Bs 


)  (A.  7) 


From  (A.5),  (A.6)  and  (A. 7),  (A.l)  becomes 


EtCPiE  Pil)  ]  —  E[(t|qB  PqL )  ]  —  2(E[r|IE  ]  E[r|ltr|IL]) 

2  P. [dM  +  ' «M  - i)(l  +  d - i /  k)  Si”(i"B^  1B-  A  2  S (M  - i)(l  - i  / k)  7  B- > 


inBf  /  Bs 

-‘x‘(M-i)(i-|kd-i|/t)sin(';B'7B-)] 


i7iBf  /  Bs 


(A.  8) 


ijiBf  /Bs 

Based  on  Central  Limit  Theorem  (CLT),  it’s  obvious  that  all  the  noise  components  and  r^are  zero 


mean  Gaussian  random  variables.  It  is  known  that,  for  zero  mean  Gaussian  variables  r|i;r|2,r|3,r|4, 

E['n1Ti2ri3ri4]  =  E[ri1ri2]E[ri3-n4]  +E[r|,ri3]E[r|2r|4]  + E[r|,r|4]E[r|2ri3]  (A.9) 

Therefore,  we  have 


f 1 9  ip  (h  ie  h  [i. )  ]  — Efriqp  (t]qE  Pql)  ] 

=  2(E[r|II.  t|ie  ]  E[r)  £p  r)iEr|]L]) 

=  2(E[r|  IP2  ]E[r|  1E2  ]  +  2E2  [x]  lpr|  IE  ]  -  E[r|  lp2  ]E[p  ffir|IL  ]  -  2E[r|  Ipr|IE  ]E[ri  IPr|IL  ]) 
-  2(E[r|  jp 2  ]E[r|  IE  2  ]  -  E[r|  ip  2  ]E[r|  IEr|  IL  ]) 

=  2E[r|IP  ](E[r)  IE  ]  E[r|  IE  rj IL  ]) 


k-l 

=  2P„2[M  +  2£(M-i)(l-i/k) 

i=l 


sin(i7rBf  /Bs) 
i7iBf  /Bs 


k-kd-l 

-][dM  +  Y  (M  —  i)(l  +  d  —  i / k) 

i=i 


sin(i7iBf  /  B s ) 
inB  f  /  B  s 


+  2£(M-i)(l-i/k) 

k-kd 


sin(i7iBf  /Bs) 
i7iBf  /Bs 


k+kd-l 

£  (M  —  i)(l  —  |kd  —  i|  /  k) 


i-[ 


sin(inBf  /  Ba ) 
inB f  /Bs 


] 


(A.  10) 
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Figure  1.  (a)  PDF  of  the  simplified  Middleton  noise  model,  (b)  Middleton  noise  sequence  with  unit 

variance 
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Figure  2.  (a)  PDF  of  generalized  Cauchy  noise  model,  (b)  Generalized  Cauchy  noise  sequence  with 

unit  variance. 
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Figure  3.  (a)  PDF  of  the  power  peaks,  (b)  PDF  of  the  low  power  values,  (c)  PDF  of  the  inter-arrival 
times  between  successive  peaks,  (d)  Sample  sequence  of  ignition  noise  with  unit  variance 
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Figure  4  (a).  PDF  of  the  amplitude  of  aggregate  UWB  signals,  (b)  Sample  sequence  of  UWB  noise 

with  unit  variance. 
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Figure  5.  GPS  DLL  cross-correlation  process 
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(a)  (b) 


Figure  6.  (a)  Early  and  late  correlation  functions,  (b)  Discriminator  function  with  1  chip  spacing 

and  2  samples  per  chip 


Figure  7.  The  frequency  spectrum  of  C/A  code  with  2  MHz  bandwidth  Butterworth  precorrelation 

filtering  (a)  before  filtered,  (b)  after  filtered 
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Figure  8.  (a)  Early  and  late  correlation  functions,  (b)  Discriminator  function  with  1  chip  spacing 

and  2  MHz  precorrelation  filtering 
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Figure  9.  (a)  Early  and  late  correlation  functions,  (b)  Discriminator  function  with  1  chip  spacing 
and  2  samples  per  chip  under  -10  dB  impulsive  noise 
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Figure  10.  The  discriminator  error  variance  through  different  bandwidth  precorrelation  filter  (a) 
under  UWB  noise,  (b)  under  Middleton’s  impulsive  noise 
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Figure  11.  The  components  of  discriminator  error  variance  through  different  bandwidth 
precorrelation  filter  (a)  under  UWB  noise,  (b)  under  Middleton’s  impulsive  noise 
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Figure  12.  The  discriminator  error  variance  with  different  sample  rate  (a)  under  UWB  noise,  (b) 

under  Middleton’s  impulsive  noise 
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Figure  13.  The  components  of  discriminator  error  variance  with  different  sample  rate  (a)  under 
UWB  noise,  (b)  under  Middleton’s  impulsive  noise 
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Figure  14.  The  components  of  discriminator  error  variance  at  different  SNR  (a)  under  UWB  noise, 

(b)  under  Middleton’s  impulsive  noise 
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Chapter  5 

Maximum  Signal-to-Noise  Ratio  GPS  Anti-Jam  Receiver 

with  Subspace  Tracking 

1.  Introduction 

Global  Positioning  System  (GPS)  is  a  tool  to  determine  position,  velocity,  and  precise  time  worldwide  by 
measuring  the  time-of-arrival  of  signals  emitted  from  satellites.  In  addition  to  its  original  military 
purpose,  GPS  has  found  a  wide  range  of  civilian  applications  such  as  navigation,  land  surveying  and 
mapping,  and  timing  and  synchronization  for  telecommunication  networks. 

For  GPS  applications,  the  main  challenges  are  the  vulnerability  of  the  GPS  receivers  to  strong 
interference  and  the  multipath  effects  on  receiver  synchronization.  GPS  employs  spread-spectrum  (SS) 
signaling,  which  provides  a  certain  degree  of  protection  against  interference.  However,  if  the  interfering 
signal’s  power  exceeds  the  30  dB  processing  gain  offered  via  the  spreading/dispreading  of  the  GPS  C/A 
signal,  the  receiver  is  unable  to  recover  the  navigation  information  conveyed  in  GPS  signal.  Therefore, 
the  design  of  GPS  receivers  must  mitigate  the  interference  and  combat  its  effect  on  the  receiver’s  ability 
to  synchronize  with  different  satellites.  Multipath,  on  the  other  hand,  is  caused  by  signal  reflections  and 
diffractions  between  the  satellite  and  the  GPS  receiver.  In  GPS,  the  desired  signal  is  the  direct  path  signal. 
All  other  signals  distort  the  desired  signal  and  lead  to  ranging  measurement  errors. 

In  this  chapter,  we  propose  an  interference  suppression  scheme  which  combines  subspace  tracking  and 
adaptive  beamforming.  Specifically,  the  received  signal  is  first  projected  into  its  noise  subspace.  The 
resulting  interference-free  signal  is  then  processed  by  a  spatial  filter,  whose  weights  are  determined  by  the 
maximum  signal-to-noise  ratio  (MSNR)  criterion.  Computer  simulations  have  shown  that  the  proposed 
method  is  effective  in  combating  strong  interference  and  enhancing  the  GPS  signal. 
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2.  Subspace  Tracking  Interference  Suppression 


2.1.  Signal  Model 

The  GPS  receiver  is  equipped  with  an  M-element  spatial  array,  as  shown  in  Figure  1.  The  waveforms 
impinging  on  the  array  are  those  of  the  GPS  signal  and  its  multipath,  interference,  and  noise.  After  down- 
conversion  and  chip-rate  sampling,  the  received  signal  vector  from  the  antenna  array  can  be  presented  in 
discrete-time  format  as 

K  L 

x(n)  =  Y,sk  ( n)ck  ( nTs -  rk  (n»a,  +  J\/  ,(/!)<!, +v(«)  (5.1) 

k= 0  1=1 

where 

T  Nyquist  sampling  interval; 

K  number  of  multipath  components; 

s  (n)  /fth  signal  component; 

ck  /fth  C/A-code  sample; 

zk  (n)  time-delay  of  the  Mi  component; 

ak  spatial  signature  of  the  /fth  satellite  multipath; 

L  number  of  interferes; 

u.  (n  )  waveform  of  the  /th  interferer; 

d/  spatial  signature  of  the  /th  interferer; 

v(«)  additive  white  Gaussian  noise  sample  vector. 

Due  to  the  weak  cross-correlation  of  the  C/A-codes,  only  one  satellite  is  considered  in  Eq.  (5.1).  The 
subscript  0  is  designated  to  the  direct-path  signal.  Let  s(«)  =  s0(n)c0(nTs  -r0(n))a0  denote  the  data 
vector  across  the  array  due  to  the  direct-path  signal.  Then,  Eq.  (5.1)  can  be  rewritten  as 
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x(//)Q  s(/;)  +  sn  (/;)  +  n(  // )  +  x(ri) 


(5.2) 


where  sm(«)  denotes  the  contributions  from  K  multipath  reflections, 

K 

Sm(«)  D  YjSk(n)Ck(nTs  ~  (5-3) 

k= 1 

L 

and  u(n)  □  y^  n;(/?)d;  is  the  compound  interference  vector. 

i= i 

2.2.  Subspace  Tracking  Based  Interference  Suppression 

Under  the  assumption  that  the  GPS  signals,  interference,  and  noise  are  independent,  the  covariance  matrix 
of  the  received  signal  becomes 

Rv*  =  E  {x(«)x"  («)}  =  Rs  +  R„  +  R,  (5.4) 

where  E{-}  represents  the  statistical  expectation,  (  )"  denotes  conjugate  transpose,  and  Rs,  Ru ,  and 

R,  are  the  covariance  matrices  of  the  GPS  signals,  the  interference,  and  the  noise,  which  are  defined, 
respectively,  as: 


RsDJE{[s(n)  +  sm(n)][s(n)  +  sm(n)f} 

(5.5) 

R0B  A'  {u(«)u,/(«)} 

(5.6) 

R  vD£{v(n)vH(fl)}  =  (r\ 

(5.7) 

where  IM  is  an  MxM  identity  matrix. 

The  subspace  tracking  based  GPS  anti-jam  receiver  is  motivated  by  the  fact  that  in  GPS,  the  desired 
GPS  signals  are  well  below  the  noise  floor  (usually  20  to  30  dB  below  the  noise  floor).  As  such,  the  total 
received  signal  power  is  dominated  by  the  jamming  signals.  In  this  case,  the  covariance  matrix  Rvv  is 
approximated  as  [1] 

R  « R  +R  (5.8) 

XX  U  V  V  ' 
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By  performing  singular  value  decomposion  (SVD)  of  Rw,  we  can  effectively  decompose  the  received 


signal  into  two  subspaces: 


ML  M 

R«  =  !>,<  ~  +a;  £  e,ef 


1=1 


□  U^jUf  +UVLVU" 


(5.9) 


where  L,  =  diag  j  2, ,  •  •  • ,  XL }  is  an  LxL  diagonal  matrix  whose  elements  are  the  L  largest  eigenvalues. 

Ui  is  an  MxL  matrix  whose  columns,  eigenvectors  e, ,  •  •  ■  eA  associated  with  L  largest  eigenvalues,  span 
the  interference  subspace.  Ev  =  cr*IM_L  contains  the  rest  M-L  eigenvalues,  which  are  assumed  to  be 
equivalent  to  av  ,  and  the  columns  of  the  MxM-L  matrix  Uy  are  the  associated  M-L  eigenvectors,  which 
span  the  noise  subspace.  Note  that  vectors  jd, ,  •  •  •  dA}  also  span  the  interference  subspace,  i.e., 

spanjej,  •••,  ei}  =  span{d1,  •••,  dz}  (5.10) 

Even  though  the  spatial  signatures  d, ,  l=\,—f  of  the  interferers  are  usually  unknown  at  the  receiver, 
the  interference  subspace  can  be  explicitly  obtained  by  first  computing  an  estimate  of  Rn  ,  then 

performing  SVD.  For  large  arrays,  however,  the  eigendecomposition  imposes  heavy  computational 
burdens  on  the  receiver.  Therefore,  such  method  may  not  be  suitable  for  real  time  processing  of  GPS 
signals.  An  alternative  approach  is  to  use  the  so  called  subspace  tracking  techniques  [2],  Such  techniques 
are  especially  suitable  for  GPS  because  the  power  of  the  interferer  is  much  stronger  than  that  of  the  GPS 
signal,  and  the  number  of  interferers  is  limited  which,  in  turn,  limits  the  dimension  of  the  interference 
subspace.  Subspace  tracking  estimates  the  interference  subspace  recursively  on  a  sample-by-sample  basis 
and,  thus,  avoids  explicit  calculation  of  the  matrix  Rn  .  For  the  anti-jam  GPS  receiver,  we  use  the 

projection  approximation  subspace  tracking  with  deflation  (PASTd)  method  proposed  in  [3].  The  PASTd 
algorithm  first  estimates  the  most  dominant  eigenvector  and  the  projection  of  the  received  data  onto  this 
eigenvector  is  then  removed  from  the  received  data.  Now  the  second  dominant  eigenvector  becomes  the 
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most  dominant  one  in  the  updated  data  vector  and  it  can  be  extracted  as  well.  Repeating  this  procedure, 
all  desired  eigencomponents  can  be  estimated  sequentially.  The  PASTd  algorithm  is  summarized  in 
Table  1. 


Table  1:  The  PASTd  Algorithm 


Usually,  the  initial  e(0)  can  be  chosen  from  an  identity  matrix.  In  the  above  algorithm,  e\  is  an 
estimate  of  the  /th  eigenvector  and  yl  is  the  corresponding  eigenvalue,  and  0  <  |3  <  1  is  the  forgetting 
factor. 

Once  the  interference  subspace  is  available,  the  noise  subspace  can  be  obtained  from  the  orthogonal 
projection  of  the  interference  subspace,  which  is  given  by 

U^=IM— U^UXr'Uf  (5.11) 

where  (□)  *  denotes  matrix  inverse.  Therefore,  columns  of  Uj1  span  the  noise  subspace.  The  projection 
of  x(«)  onto  Uj1  yields 

y(«)  =  U^x(h)  =  [s(«)  +  sm(»)]  +  Uf  v(w)  (5.12) 

which  only  contains  contributions  from  the  GPS  components  and  noise. 
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3.  MSNR  Beamformer 


From  the  above  discussion,  we  know  that  by  projecting  the  received  data  onto  the  noise  subspace,  the 
interfering  signals  are  completely  suppressed.  After  the  suppression  of  the  interference,  the  GPS  signal  is 
still  far  below  the  noise  floor.  In  order  to  synchronize  the  receiver  with  the  satellite,  which  is  usually 
achieved  by  cross-correlating  the  received  data  with  a  locally  generated  C/A-code  and  identifying  the 
maximum  value,  the  GPS  signal  must  be  enhanced.  To  this  end,  we  design  a  filter  such  that  the  output  of 
the  filter  achieves  the  maximum  signal-to-noise  ratio  (MSNR). 

Let  w  be  the  Mx  1  weight  vector.  Then,  the  output  of  the  filter  is  given  by 

z(n )  =  w^y  (n) 

(5.13) 

=  [s(n)  +  sm  («)]  +  w//UILv(«) 


and  the  filter  w  is  determined  from 


wmsnr  —  tnax  - 


E\  w"U^[s(»)  +  sm(n)] 


£  w  Uj  \(n) 


(5.14) 


=  max 


ff  U,R,Ui  w 
w"UfUf"w 


Equation  (5.14)  indicates  that  the  determination  of  beamformer  w  requires  the  power  of  the  GPS 
signal  (contains  both  the  power  from  the  direct-path  signal  and  the  contributions  from  its  multipath 
components)  at  the  receiver.  Generally,  the  calculation  of  the  GPS  signal  power  requires  some  a  priori 
knowledge  of  the  satellite.  For  example,  as  indicated  in  [4],  if  the  location  of  the  satellite  is  known,  then 
the  GPS  signal  power  can  be  computed  as  follows.  Assume  that  the  satellite  is  located  at  the  angle  ( 6 \  y/). 
For  GPS,  the  temporal  autocorrelation  of  the  transmitted  signal,  which  is  essentially  the  autocorrelation 
function  of  the  Gold  code  denoted  as  R(  (r)  is  known.  The  MxM  matrix  R  is  calculated  in  the  absence 


of  interference  and  noise.  Let  77(77?)  denote  the  phase  shift  for  the  GPS  satellite  at  the  angle  (0,\|/)  to  the 
777th  antenna  and  let  r  □  f  J  .  Then,  R  is  given  by  R  =  (0)rr7/ . 
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On  the  other  hand,  if  the  satellite  location  information  is  not  available  at  the  receiver,  we  take  an 
alternative  approach  to  solve  the  above  maximization  problem.  From  Eq.  (5.12)  we  know  that  y (n)  only 


contains  contributions  from  the  GPS  signals  and  noise.  Due  to  the  weakness  of  the  GPS  signals,  the 
output  of  the  projection  is  dominated  by  noise.  Note  that 


E 

jw"y(«) 

I  _  w"UltR,,UlL"w 

£{ 

w//UILv(«) 

21 

J 

|  <72Wi/UjLUjU/W 

w"U^RsU^w 
a;w"UfU^w  + 


which  shows  that  the  beamformer  w  that  maximizes 


w^U^R.U^w 


also  maximizes 


wtfUfR,,U^w 

o-v2w"UfUf"w 


Therefore,  the  filter  w  can  be  found  by  solving 


WMSNR 


=  max 

w 


o-;w"u^u^w 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


Thus,  the  optimum  w  is  the  eigenvector  corresponding  to  the  dominant  eigenvalue  of  the  following 
generalized  eigenvalue  problem: 


UXXWw  =  //U^U]L//w  (5.19) 

where  //  denotes  the  dominant  eigenvalue,  which  is  also  the  maximum  SNR. 

In  practice,  Rtx  is  replaced  by  its  sample  estimate,  which  can  also  be  computed  recursively  [5]. 
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4.  Simulations 


A  linear  uniform  array  consisting  of  M=1  sensors  with  half-wavelength  spacing  is  used  in  the  simulation. 

In  the  first  simulation,  we  investigate  the  convergence  of  the  PASTd  algorithm.  In  this  experiment, 
there  is  one  jammer  added  to  the  received  signal  and  we  set  the  signal-to-noise  ratio  (SNR)  to  -10  dB  and 
signal-to-interference  ratio  (SIR)  to  -30  dB.  This  lead  to  the  interference-to-noise  ratio  (INR)  of  20  dB. 
Figure  2  shows  that  after  250  samples,  the  eigenvalue  converges  to  the  true  INR  of  20  dB. 

In  the  next  experiment,  we  examine  the  interference  suppression  performance  of  the  proposed  GPS 

o  o  o 

receiver.  There  are  two  jammers  located  at  30  and  60  ,  and  the  satellite  is  at  10  .  Figure  3  shows  that  the 
receiver  can  successfully  generate  high  gain  toward  the  satellite  direction,  while  placing  deep  nulls  at  the 
jammer  locations.  By  using  the  proposed  GPS  receiver,  the  desired  GPS  signal  is  enhanced,  whereas  the 
jammers  are  suppressed. 

Finally,  we  investigate  the  proposed  receiver’s  synchronization  capability.  Generally,  the 
synchronization  can  be  achieved  by  cross-correlating  the  received  signal  with  the  locally  generated  C/A- 
code  [6].  When  the  receiver  synchronizes  with  the  satellite,  there  is  a  maximum  correlation.  In  the 
simulation,  we  considered  three  different  scenarios  to  illustrate  the  receiver’s  performance.  We  set  SNR= 
-30  dB  and  SIR=- 40  dB.  In  the  first  case,  the  received  signal  is  directly  correlated  with  the  C/A-code  and 
the  resulting  normalized  cross-correlation  is  shown  in  Figure  4(a).  It  is  noted  that  without  any  processing, 
the  synchronization  fails.  If  only  the  interference  suppression  is  applied,  the  receiver  is  able  to 
synchronize  with  the  satellite  after  the  cross-correlation,  but  the  noise  contribution  remains  significant,  as 
shown  in  4(b).  With  the  proposed  receiver,  however,  the  noise  can  be  drastically  reduced.  This  is  shown 
in  4(c). 
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5.  Conclusions 


In  this  chapter,  we  have  considered  the  problem  of  interference  cancellation  in  GPS.  Specifically,  a  GPS 
receiver  combining  the  subspace  interference  suppression  and  MSNR  beamforming  is  proposed.  Through 
computer  simulations,  we  have  shown  that  the  proposed  receiver  is  capable  of  providing  high  gains  for 
the  desired  GPS  signal  while  suppressing  the  strong  interferers. 
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Figure  1.  Block  diagram  of  the  proposed  GPS  receiver. 


Figure  2.  Convergence  of  the  eigenvalue. 
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Figure  3.  Beampattern  of  the  proposed  GPS  receiver. 
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